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Abstract 

We present a derivation of a Boltzmann equation for the QCD plasma, starting 
from the quantum field equations. The derivation is based on a gauge covariant 
gradient expansion which takes consistently into account all possible dependences 
on the gauge coupling assumed to be small. We point out a limitation of the gradient 
expansion arising when the range of the interactions becomes comparable with that 
of the space-time inhomogeneities of the system. The method is first applied to the 
case of scalar electrodynamics, and then to the description of long wavelength colour 
fluctuations in the QCD plasma. In the latter case, we recover Bodecker's effective 
theory and its recent reformulation by Arnold, Son and Yaffe. We discuss interesting 
cancellations among various collision terms, which occur in the calculation of most 
transport coefficients, but not in that of the quasiparticle lifetime, or in that of the 
relaxation time of colour excitations. 



CERN-TH/99-71 
March 1999 

Published in Nuclear Physics B 557 (1999) 183-236. 



a E-mail : blaizot @spht . saclay. cea. fr 

6 Laboratoire de la Direction des Sciences de la Matiere du Commissariat a l'Energie Atomique 
c E-mail : edmond . iancu@ccrn . ch 



1 Introduction 



At high temperature, non-Abelian gauge theories describe weakly coupled plasmas whose 
constituents, quarks and gluons for hot QCD, have typical momenta k ~ T, where T is the 
temperature |2|, |3[| . The plasma particles may take part in collective excitations which 
develop typically on a space-time scale A ~ 1/ gT much larger than the mean interparticle 
distance f ~ 1/T (g is the gauge coupling, assumed to be small). Such collective exci- 
tations can be described in terms of mean fields carrying appropriate quantum numbers 
and coupled to the plasma particles. In this long- wavelength limit, the plasma particles 
obey simple, collisionless, kinetic equations which can be viewed as the generalization of 
the Vlasov equation of ordinary plasmas ||, [5], [| . 

In this paper, we shall be interested in specific collective excitations involving colour 
fluctuations on larger wavelengths, A ~ 1 / g 2 T. In this situation, the effects of the colli- 
sions among the plasma particles become as important as those of the mean fields. The 
kinetic equations obeyed by the plasma particles must therefore be generalized so as to 
include the collisions terms. This is what we shall do here. The Boltzmann equation 
that we shall obtain (see eqs. ( |1.1| ) — ( |1.3| ) below) turns out to be identical to the one 
proposed recently by Arnold, Son and Yaffe 0, and yields, in leading logarithmic accu- 
racy, Bodeker's effective theory for the soft (p ~ g 2 T) fields §. The derivation presented 
here, starting from the quantum field equations, clarifies the nature of the approxima- 
tions involved, and thus fixes its range of applicability. Furthermore, it also provides 
some justification for numerous previous works using ad hoc transport equations inspired 
by classical transport theory || |9], [It], [11], [12], [L3|, [L4|, [L5], |IEfl . It should be emphasized 
that transport equations with a similar colour structure have been also proposed in Refs. 
17 , [To[| , and that some of the technics that we shall be using have been used already by 



many authors [Tj| [T^, |2(| ^T], 2"2"|| . However, in most of these works, the organizing role of 
the various dynamical scales which appear in hot QCD was not recognized, which led to 
unecessarily complicated, and sometimes inconsistent, equations. 

In fact, "kinetic theory" in the way we use it here could be regarded as a powerful 
tool for constructing effective theories for the soft modes of the plasma. These soft degrees 
of freeedom are represented by mean fields, while the hard ones, which are "integrated out" 
using perturbation theory, survive as induced sources for these mean fields. The resulting 
effective theory can then be studied non-perturbatively, e.g., as a classical theory on a 
lattice: recently, this strategy has received much attention in connection with studies of 
baryon number violation in the high temperature electroweak theory |23|, |53, [25], |26], [| 0, 
p7| , |T6|| . Let us also recall that this method has been first demonstrated for the collective 
dynamics at the scale gT, where we have shown [f|, || that simple, collisionless, kinetic 
equations resum an infinite number of one-loop diagrams with soft external lines and hard 
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loop momenta, the so-called "hard thermal loops" pB, 25 1. 

Let us now summarize the main equations to be obtained below. The collective, 
longwavelength colour fluctuations of the hard transverse gluons are described by a density 
matrix N(k,x) which, to the order of interest, can be written in the form: 

iV ab (k, x) = N(e k )S ab - gW ab (x, v) (diV/de fc ), (1.1) 

where N(ek) = l/(e* £k — 1) is the Bose-Einstein thermal distribution (with = |k|), and 
the function W(x, v), which parametrizes the off-equilibrium deviation, is a colour matrix 
in the adjoint representation, W(x,v) = W a (x, v)T a , which depends upon the velocity 
v = k/^fc (a unit vector), but not upon the magnitude k = |k| of the momentum. The 
functions W a (x, v) satisfy the following transport equation: 

(v ■ D x )°»W b (x,v) = v ■ E B (x) - 7 {w a (x, v) - ^ ^ } ■ (1.2) 

In the left hand side of this equation, v ■ D x is a gauge- covariant drift operator (with 
= (l,v) and = <9 M + igA 11 ), while in the right hand side we recognize a mean 
field term (v • E a (a;), with E a the chromoelectric field) and a collision term. The latter is 
proportional to the quasiparticle damping rate 7, which appears to set the scale for the 
colour relaxation time: t co \ ~ I/7 ~ l/(g 2 T ln(l/g)) (see below). The other notations 
above are as follows: the angular brackets in the collision term denote angular average 
over the directions of the unit vector v' (as in eq. (p. .4|) below), and the quantity $(v ■ v') 
is given by: 

$(v • v') = (2vr) 2 J 5(q - q ■ v)%o ~ Q ' V)|*A( g ) + (v t ■ v' t ) *V t (q)\* . (1.3) 

where *T>i(q) and *T> t (q) denote the resummed gluon propagators in the electric and 
the magnetic channels, respectively |l], 0- Up to a normalization, $(v • v') is the total 
interaction rate for two hard particles with momenta k and p (and velocities v = k and 
v' = p) in the (resummed) Born approximation, as illustrated in Fig. [1] (vt and are the 
transverse projections of the velocities with respect to the momentum q of the exchanged 
gluon: e.g., v\ = (<5 U — q l q>)vi). The damping rate 7 is obtained from $(v • v') as follows: 

3 4 iV c 2 T 3 /-dO' . „ g 2 N c T 
7 = y —^- J — $(v • vO ^ y -^- Hl/g\ (1.4) 

and is of 0(g 2 T) in spite of the explicit factor g 4 in front of the above integral. This 
is because the quasiparticle damping is dominated by soft momentum transfers q ^ gT, 
which gives an enhancement factor ~ 1 / g 2 after the resummation of the screening effects 
at the scale gT @, |3T| |3T], g2|. 
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p'=p+q 



k k'=k-q 

Figure 1: Elastic scattering in the (resummed) Born approximation. The continuous lines 
refer to hard gluons (these are off-equilibrium propagators), while the wavy line is the 
soft gluon exchanged in the collision. The blob stands for HTL resummation. 



Actually, in the present approximation, 7 is even logarithmically infrared divergent, 
due to the unscreened static (go 0) magnetic interactions. (In writing the right hand 
side of eq. ( |L4 ) we have assumed an infrared cutoff ~ g 2 T, as it is usually done in the 
literature [ ST] , ID , [T2|].) To logarithmic accuracy, that is, by preserving only the singular 
piece of the magnetic scattering element in eq. (|1.3| ) (see eq. ( |3.123|) below), eqs. (|1.2j )- 
( |1.3|) generate Bodeker's effective theory for the soft modes ||. 

Within the same accuracy, eq. ( |1.2| ) can be solved to get the so-called colour con- 
ductivity |LC], [L2], || 0. The induced colour current is expressed in terms of W(x, v) 
as: 



j a (x) = m 2 D (vW a {x^)), 



(1.5) 



with the Debye mass m? D = g 2 N c T 2 /3. For constant colour electric fields, we get from 
eq. (0): 

1 
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so that 
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with the colour conductivity a c = m^/37 ~ T j \n(l/g). 

The next section of the paper contains a derivation of the Boltzmann equation for 
scalar electrodynamics (SQED). There are several reasons for this. First, to our knowl- 
edge, this is the first consistent derivation of a Boltzmann equation for gauge theories, 
starting form the quantum field equations. Second, it serves as a preparation for the more 



3 



involved non Abelian case of QCD which is presented in the following section. Finally, 
and this is the most important, it will reveal interesting compensations which occur in 
Abelian, but not in non- Abelian gauge theories. Thus, in SQED, most transport phe- 
nomena are dominated by large angle scattering, so that the typical relaxation times are 
r tr ~ l/(e 4 Tln(l/e)), where e is the electric charge || [15|]. This is to be contrasted 
with the quasiparticle lifetimes which are limited by small angle scatterings and are of 
order r ~ 1/ (g 2 T\n(l/g)) (in both QED and QCD). The same cancellations occur in most 
cases for QCD as well (thus yielding, e.g., a viscous relaxation time r„ sc ~ l/(g 4 Tln(l/g)) 
H [14]]), except for the the relaxation of colour excitations which remains dominated by 
very soft gluon exchanges [JTO], [12], §|, As a result, the colour relaxation time turns 
out to be of the same order as the quasiparticle lifetime, as is evident in eq. (|1.2|) . This 
yields a colour conductivity o c ~ T/ln(l/g), to be contrasted with the usual, electric 
conductivity^ a el ~ T/(e 2 ln(l/e)) p[. 

The main part of the paper is section 3, which contains the derivation of the Boltz- 
mann equation for the QCD plasma, and a discussion of the approximations which are 
needed in this derivation. These involve a gradient expansion of the Dyson- Schwinger 
equations, supplemented by a perturbative evaluation of the collision terms and a lin- 
earization with respect to the off-equilibrium fluctuations. All these approximations are 
commonly used in deriving kinetic equations from quantum field theories, and they are 
generally seen as independent approximations (to some extent they remain so in the case 
of SQED discussed in Sec. 2 below). However, in order to fulfill the constraints imposed 
by a non Abelian gauge symmetry, we shall see that it is convenient to control all these 
approximations by the same small parameter, namely the gauge coupling g. Thus, for 
instance, the amplitudes of the mean fields will be restricted so that \A%\ ~ gT : this 
guarantees that the two terms in the soft covariant derivative D x = d x + igA^ are of the 
same order in g, dx ~ gA ~ g 2 T, so that Dx = 0(g 2 T) can be preserved consistently 
in the expansion. A further difficulty that we shall have to face is related to the poor 
convergence of the gradient expansion when the range of the interactions becomes com- 
parable to the scale of the system inhomogeneities. As we shall see this will be the main 
limitation of the accuracy of the collision term. 

For completeness, we present in section 4 some diagrammatic interpretation of the 
Boltzmann equation. (Previously, the connection between Feynman graphs and the Boltz- 



mann equation has been explored in detail only for a scalar field theory, in Refs. [|33|, |34 
The section 5 summarizes the conclusions. 



d We mean here, of course, the electric conductivity in a QED or OCD plasma, that is, in a gauge theory 
without electrically charged vector bosons. The situation would be different in the electroweak theory 
where, in the high-temperature, symmetric, phase, the electric charge can be efficiently randomized via 



small angle scatterings mediated by the W -bosons 15 
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2 Scalar QED 



In this section, we briefly summarize the general formalism which allows one to construct 
kinetic equations from the Dyson-Schwinger equations obeyed by the non-equilibrium 



Green's functions [[35], |36], |37|, [39], f4(| |4l|. In order to bring out the essential aspects 
of the formalism while avoiding the complications specific to non Abelian gauge theories, 
we shall consider here scalar electrodynamics (SQED), with Lagrangian: 

c = (D^iDyy-^F^, (2.i) 

where is a complex scalar field, is the photon field, = <9 M + ieA^ is a covariant 
derivative, and F^ u the field strength tensor, F^ v = d^A u — d v A^. 

The systems that we consider are assumed to be initially in thermal equilibrium, 
and described by the density operator p = -|exp{ — j3H} where H is the Hamiltonian 
corresponding to (|2.1|), and Z is the partition function. At some time t , a time-dependent 
external perturbation (an electromagnetic current j M (x)) starts acting on the system, so 
that the Hamiltonian becomes: 

H 3 (t) = H + J d 3 xj(t,x)-A(x), (2.2) 

where j ■ A = A^. The density operator at time t is given by: 

p(t) = U(t,t )pU(t ,t), (2.3) 

where U(t,to), the evolution operator, satisfies: 

id t U{t,t ) = HjtyUfato), U{t ,t ) = 1. (2.4) 

In the presence of the perturbation, the gauge field A^ develops an expectation value: 

C p -f3H ~) 

Tr(p(t)^) = Tr(pA»(tj) = Tr| — A"(f)J = (A»{t)) , (2.5) 

with 

A»(t) = U- l (t,t )A^U(t,t Q ) = U(t ,t)A»U(t,t ), (2.6) 

where we have used the fact that U(t ,t) = £/ _1 (t, t ). More generally, we shall be 
interested in various n-point functions, and in particular in 2-point functions for which 
we shall derive equations of motion in the next subsection. For instance, the time-ordered 
2-point function of the charged scalar field is given by: 

G(t 1 ,t 2 ) = (T0(t 1 )0 t (t2))=Tr|^T0(t 1 )0t(t 2 )l 

= B{t l -t 2 )G > {t l ,t 2 ) + e{t 2 -t 1 )G < {t l ,t 2 ) (2.7) 
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where <j>(t) = U(t Q ,t) (fiU(t,t ), and the functions G > and G < are defined by: 

G>(x,y) = (<P(x)^(y)), G<{x,y) = (tf{y)<f>(x)). (2.8) 



(To lighten the notation, the spatial coordinates have not been indicated in eq. (|2.7| ).) 
The functions G > and G < can be used to construct the retarded (Gr) and advanced {Ga) 
propagators, which will also be needed: 

G R (x,y) = i9(x - y )[G > (x,y) - G < (x,y)\ , 

G A {x,y) = -iflfoo-Zo) 3/) - G<{x,y)]. (2.9) 

Similar definitions hold for the photon 2-point functions. In the case of the photon, 
we shall decompose the gauge field into its average value for which we shall reserve the 
notation A^x) (i.e., in the following we identify (A^) — > A^), and a fluctuating part a^x) 
with (a M ) = 0. The time ordered photon propagator is then given by: 

V^{x,y) = (Ta^(x)a u (y)) . (2.10) 

The 2-point functions introduced above satisfy boundary conditions which follow from 
their definitions (cf. eq. ( |2.7| )). For instance: 

G<(t ,z) = G > (t -if3,z), (2.11) 

and similarly for the photon 2-point functions and for the various self-energies to be 
introduced later. Furthermore, these functions have hermiticity properties which will 
be useful below. Specifically, all the "bigger" (>) and "lesser" (<) 2-point functions 
are hermitian: for instance, (G > (x, y))* = G > (y,x) and (V> u (x,y))* = V>^(y,x). This, 
together with the definitions ( p.9|) , imply (Gr(x, y))* = GA(y,x) and (Vr(x, y))* = 
V A ^(y,x), together with similar properties for the various self-energies. 

In thermal equilibrium, the system is homogeneous (e.g., Gf g (x,y) = Gf q (x — y)), 
and it is convenient to go to momentum space. Then, the boundary condition (|2.11| ) 
translates into the so-called KMS condition : 

G> q (k) = e^G<(k), (2.12) 

which implies the following structure for the equilibrium 2-point functions: 

G> (k) = P (k) [l + N(k )} , G<(k) = p(k) N(k ) , (2.13) 

with N(ko) = l/(e@ ko — 1) and the spectral density p(k) = G> (k) —GfJk). In particular, 
for free, massless particles: 

G<(k) = po(k)N(k ), G>{k) = po(k)[l + N{k )\, (2.14) 

with p (k) = 2Tre(k )5(k 2 ). 
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Figure 2: Complex-time contour for the evaluation of the thermal expectation values: 
C = C+UC-U c . 

2.1 Equations of motion for Green's functions 



In order to obtain the equations of motion for the 2-point functions, it is convenient 
to extend their definition by allowing the time variables to take complex values. More 
specifically, we introduce, in the complex time plane, the oriented contour depicted in 
Fig. |[ This may be seen as the juxtaposition of three pieces: C = C + UC_U C . We call 
z the (complex) time variable along the contour, and reserve the notation t for real times. 
On C + , z — t takes all the real values between to to tf. On C_, we set z = t — irj (r/ — > 0+) 
and t runs backward from tf to t . Finally, on C , z = t — ir, with < r < (3. We define 
a contour ^-function 9 C : @c( z i, z 2) = 1 if 2i is further than z 2 along the contour (we then 
write Z\ y z 2 ), while 9 c (zi,z 2 ) = if the opposite situation holds (z\ -< z 2 ). We can 
formalize this by introducing a real parameter u which is continously increasing along the 
contour; then, the contour C is specified by a function z(u), and 0c(zi,z 2 ) = 9{u\ — u 2 ). 
We shall later need also a contour delta function, which we define by: 



Sc{zi,z 2 ) 




5(ui - u 2 ). 



(2.15) 



The definition of the propagators is then extended in a natural way. For instance, the 
contour-ordered propagator of the scalar field becomes: 



G{ Zl ,z 2 ) = (T c 0(^)0t(^ 2 )) =Tr 



-0H 



■T c <j>(zi) 



^2 



(2.16) 



where Tc orders the operators on its right, from right to left in increasing order of the 
arguments Uj. For time arguments t ly t 2 on C + , the contour propagator ( |2.16| ) reduces to 



7 



the time-ordered propagator fl2.7p. For t\ G C_ and t 2 G C + , we have G(ti — ir},tz) = 
G > (ti,t 2 ), while for t\ G C + and t 2 G C_, we have G(ti,t 2 — 277) = G r< (ti,t 2 ). 

With these definitions in hand, most of the formal manipulations familiar in equi- 
librium field theory can be extended to the case of non equilibrium. This is convenient 
for the derivation of the equations of motion for the n-point functions to which we now 
turn. 

The mean field equation is: 

^(x)=/(x)+jL(x), (2.17) 

with the induced current 

IM?) = -ie((D^(x)) 1 <p(x) - <\>\x) (D»<f>(x))). (2.18) 

In this expression, = d^+ig{A^ + a^). However, in line with the approximations below, 
we can ignore the contribution of the quantum field a M in the expression of the induced 
current. This amounts to neglect the contribution of a connected 3-point function. We 
can then write: 

1U?) = - ( D $) G<(x,y)\ x=y , (2.19) 

where now, and for the rest of this section, = + igA^, with the average gauge 
potential. 

In order to calculate the induced current, we need the 2-point function G < (x,y). 
An equation of motion for this function can be obtained from the equation of motion for 
the time ordered propagator: 

-D 2 x G(x,y) - i f d 4 zS(x,z)G(z,y) = i5 c (x,y), (2.20) 

where 5c(xo,Vo) is the contour delta function, and E the scalar self-energy. The latter 
admits the following decomposition, similar to that of G, eq. (|2.7p: 

E(x,y) = -iE(x)8 c {x,y) + 9 c {x , y )S > (x, y) + 0c[yo, x )S < (x, y). (2.21) 

We have separated out a possible singular piece S (e.g., the standard tadpole diagram 
which generates a temperature-dependent mass correction [f|2|]). The non-singular com- 
ponents S > and S < obey a boundary condition similar to ( 2.11] ). In particular, in equi- 
librium, £>(Jfe) =e l3ko E<(k). 

The equations of motion in real-time for the mean field and the 2-point functions 
are obtained by letting the external time variables xq and yo take values on the real-time 
pieces of this contour, C + and C_. For xq G C + , the mean field equation is formally the 
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same as in eq. ( |2.17| ). Consider now eq. ( |2.20| ): by choosing x G C + and yo G C_, and 
by using the decompositions (|2.7|) and ( [2.2 1| ), we obtain, after some manipulations, an 
equation for G < (x, y): 



(Dl + t(x))G<(x,y) = - f d 4 z[i: R (x,z)G < (z,y) + £<(x, z) G A (z, y) 



together with a similar equation where the differential operator is acting on y: 
((Dl) 2 + t(y))G<(x,y) = - [ d 4 z[G<(x, z) E A (z, y) + G R (x, z) S < (^, y) 



(2.22) 



(2.23) 



In these equations, D 2 = D^D^, = <9 M — ieA^, and we have used the definitions (|2.9|) 
for the retarded and advanced Green's functions, together with similar definitions for E# 
and Ha- One can also obtain an equation satisfied by G R (x,y): 

(D 2 x + ±(x))G R (x,y) + J d i zZ R (x,z)G R (z,y) = 6®(x-y). (2.24) 

Note that, while the Green's functions G > and the G K and the corresponding self-energies 
are coupled by eqs. (|2.22|) -( p.22| ), the retarded Green's function G R is determined by the 
retarded self-energy alone. 

The above equations must be supplemented with some approximation scheme in 
which, for instance, the self-energy S is expressed in terms of the propagator G. Below, 
we shall use perturbation theory for this purpose. We shall refer to the above equations 
as the Kadanoff-Baym equations. They were first obtained in the framework of non- 
relativistic many-body theory Note that, in these equations, any explicit reference 
to the initial conditions and to the KMS condition has disappeared. These only enter as 
boundary conditions to be satisfied by the various Green's functions in the remote past. 
The same set of equations has been derived by Keldysh |38) to describe non-equilibrium 
evolutions of quantum systems (see also |36], [37], |40[|). 



2.2 Gauge covariant Wigner transforms 

For slowly varying off-equilibrium perturbations, the Kadanoff-Baym equations can be 
transformed into kinetic equations, as we now explain. In thermal equilibrium, the system 
is homogenous, and the two-point functions depend only on the relative coordinates s M = 
x m _ thermal particles have typical energies and momenta k ~ T. It follows 

that the 2-point functions are peaked around s M = 0, their range of variation being fixed 
by the thermal wavelength At = ~ 1/T. In what follows, we shall be interested in 
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off-equilibrium deviations which are slowly varying in space and time, over a typical scale 
A > A T . 

In order to take advantage of the assumed separation of scales between hard de- 
grees of freedom (the plasma particles), and the soft degrees of freedom (the collective 
excitations at scale A ^> At), it is convenient to introduce relative and central coordinates, 

s^ = x^~ y", X» = jZ- , (2.25) 

and to use the Wigner transforms of the 2-point functions. These are defined as Fourier 
transforms with respect to the relative coordinates s M . For instance, the Wigner transform 
of G < (x, y) is: 

G < (k,X) = Jd 4 se ik - s G < (x+'^,X-^), (2.26) 

and we shall use similar definitions for the other 2-point functions. Note that we shall 
use the same symbols for the 2-point functions and their Wigner transforms, considering 
that the different functions can be recognized from their arguments. 

The hermiticity properties of the 2-point functions discussed after eq. ( |2.11| ) imply 
similar properties for the corresponding Wigner functions. For instance, from (G > (x, y))* — 
G > (y,x) we deduce that G < (k, X) is a real function, (G < (k, X))* = G < (k,X), and simi- 
larly for G > (k, X). Also, (G^k, X))* = G^ik, X). Similar properties hold for the photon 
2-point functions and for the various self-energies. 

In gauge theories, the physical interpretation of the Wigner functions as phase 
space densities is complicated by the lack of gauge covariance of the 2-point functions. 
To remedy this, we shall define new, gauge invariant, functions, whose construction may 
be motivated by considering the conserved electromagnetic current: 

f{x) =ie(D%- (£^) f ) G<(x,y)\ x=y , (2.27) 

where G < (x,y) = (ft (y)</>(x)) is not gauge invariant. It is easy to define a corresponding 
gauge invariant function by multiplying it by a parallel transporter, or 'Wilson line" , 

U(x,y)=e- te ^ dz ^ (z \ (2.28) 

where the path 7 joining y to x is a priori arbitrary. Thus, for instance, 

G<(x,y) = (4>\y)U(y,x)4>(x)) = U(y,x)G < (x,y) (2.29) 

is manifestly gauge invariant. The conserved current may then be expressed in terms of 
this gauge invariant function: 



p{x) = (d% - d%) G<(x,y) 



x=y 



2id^G<(s,X) s=Q . (2.30) 
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To see this, note that G < (x,y) = U(x,y)G < (x,y), and that Dg U(x,y)\ x=y = 0. Note 
also that the expression for is independent of the path joining x and y, since only 

an infinitesimal path is needed. 

For definitness, we shall in fact choose 7 to be the straight line joining x and y. 
This choice is physically motivated since, as we shall see later, the hard particles preserve 
straight line trajectories in the presence of the soft mean fields (at least, to leading order 
in e). Moreover, as shown in Refs. |[8|, [L9|, ^T], such a path allows one to interpret the 
covariantization procedure as the replacement of the canonical momentum by the kinetic 
one (see eq. ( |2.32|) below). This being said, most of our results below will be independent 



of the exact form of 7 (see, however, the discussion after eq. ( |2.41| )). Indeed, we shall 
mostly need the parallel transporter U(x, y) in situations where the end points x and y 
are close to each other (\s\ ^ 1/T, with s = x — y), so that the variation of the field can 
be neglected along the path. This is a good approximation provided 7 never goes too far 
away from x and y, that is, provided \z — X\ = 0(1/T) (with X = (x + y)/2) for any 
point z on 7. For any such a path we can write: 

U{x,y) « e - ies ' A W, (2.31) 

up to terms which involve, at least, one soft derivative dxA^ (and which do depend upon 
the path). 

Starting from G < (x, y) , we construct the gauge invariant Wigner function: 

4„ Ak-s tt I v s v 1 s \ n< ' •"■ s -- s 



G<(k,X) = J d'se^U[X--,X + -j G<(X + -,X-- 
w Jd i se* s < k+eAix »G < (x+ S -,X- S -^) 

= G < (p = k + eA(X),X). (2.32) 

This formula shows that, as alluded to before, the gauge invariant Wigner function may 
be obtained from the ordinary one by the simple replacement of the canonical momentum 
p M by the kinetic momentum k^ = — eA^(X). Returning to the current, we see that 
it takes the form: 

f{x) = 2eJ-0^(p»-eA»(X))G<(p,X) 

= 2eJ-0^k»G<(k,X). (2.33) 



These two expressions for the current may be seen as the analogs of eqs. fl2.27| ) and ( |2.30| ) . 
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2.3 Gradient expansion and kinetic equations 



For slowly varying disturbances, taking place over a scale A ^> Ay, we expect the s M 
dependence of the 2-point functions to be close to that in equilibrium. Thus, typically, 
k ~ d s ~ T, while 9x~l/A<T. The general equations of motion written down in Sec. 
2.1 can then be simplified with the help of a gradient expansion, using k and X as most 
convenient variables. 

The starting point of the gradient expansion is the equation obtained by taking the 
difference of the Kadanoff-Baym equations ( |2.22|) and (|2.23|) . For further reference, we 
shall call it the difference equation. We then define: 

E(x,y)=I%- (Dj) 2 , (2.34) 

where: 

D 2 X = dl + 2ieA(x) ■ d x + ie(d ■ A{x)) -e 2 A 2 (x), 
(Dj) 2 = d 2 y -2ieA(y)-d y -ie(d-A(y))-e 2 A 2 (y). (2.35) 

By replacing the coordinates x M and y^ by s M and X^ (cf. eq. ( |2.25| )), and rewriting the 
derivatives as: 

d x = d s + ±d x , d y = -d s + \d x d 2 x -d 2 y = 2d s -d x , (2.36) 

we perform a gradient expansion in E, with d s ~ T and 9jf~l/A<CT, and preserve all 
the terms involving at most one soft derivative d x . For instance, 

A^(X + 8 /2) w A M pO + (1/2) (s • a x )A M (X). 

A straightforward calculation yields then: 

E{s, X)n2d s -d x + 2ieA„{X)d x + 2ie{s ■ d x A^ + 2ie{d x ■ A) - e 2 (s ■ d x A 2 ) + • • ■ , 

(2.37) 

where the dots stand for terms which involve at least two soft derivatives d x . 

Before taking the Wigner transform, we make the difference equation covariant by 
multiplying both sides by the parallel transport U(y,x) (cf. eq. ( |2.29 )). For the left hand 
side, we use the expansion of E(s,X), together with eq. (ggg ) to obtain: 

U(y,x) (D 2 X - (Dl) 2 ) (u(x,y)G<(x,y)) « 2(d s ■ d x + ies^(X)d^)d<(s, X). (2.38) 

The right hand side of the difference equation involves convolutions of the form: 

C(x,y) = Jd 4 zi;(x,z)G(z,y). (2.39) 
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Upon multiplication by U(y,x), this becomes the gauge invariant quantity: 

C{x,y) = U(y,x) f d 4 z U(x, z)±(x, z)U(z,y)G(z,y) 



= f ^zP{x,y,z)t{x,z)G{z,y), (2.40) 

where we have set C(x, y) = U (y, x)C(x, y), and P(x, y, z) denotes the following plaquette: 

P(x, y, z) = U(y, x)U(x, z)U(z, y). (2.41) 

In line with the approximations in eq. (|2.38| ), we need the gradient expansion of eqs. ( |2.40| ) 
and (|2.41|) up to terms involving one soft derivative of the background field. In each of 
the parallel transporters, we choose the path 7 to be the straight line (cf. the discussion 
before eq. ( |2.31| )). Then, the plaquette ( [2.41| ) can be easily expanded around the point 
X = (x + y)/2 to yield: 

P(x,y,z) «exp{-| S ^(X)^} « 1 - ^F, U (X)5\ (2.42) 

where s = x — y, X = (x + y)/2, 5 = 2(z — X). 

We are now in position to take the Wigner transform. The only delicate step 
concerns the transformation of C (x, y) , which is given by: 

C(k,X) » t(k,X)G(k,X) + % - {£, 6} pB - l -eF, u (X) {ST h 6) + ... , (2.43) 

where {A, B}p_b. denotes a Poisson bracket: 

{A, B} pB = d k A ■ d x B - d x A ■ d k B . (2.44) 

The third term in the r.h.s. of eq. ( [2.43 ), involving F^ U (X), comes from the plaquette 
([2.42 ) and is therefore sensitive to the choice of 7. To the accuracy where this term is 
important, we expect the gauge-invariant Wigner functions X) and G(k,X) to be 
path-dependent as well. However, such path-dependent terms will disappear in the final 
form of the Boltzmann equation that we shall obtain (see eq. (|2.56[) below). 

By using eqs. ( ^.38[ ) and (2.43), the difference equation finally becomes: 



2{k-d x -ek-F-d k )G < + (d x t)d^G < -{Re± R ,G < } pB -{± < ,ReG R } 



P.B. 



+ eF^ {(d k u ±<)(d k u ReG R ) + (d k Re± R )(d k u G<' 



= -(g > S < -S > G < ) (2.45) 

In deriving the equation above, we have used the following relations: 

G R (k, X) - G A (k, X) = i(G>(k, X) - G < (k, X)) = tp{k, X), 

t R {k,X)-t A {k,X) = i(± > (k,X)-± < (k,X)) = -iF(k,X), (2.46) 
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which follow, e.g., from the definitions ( |2.9| ) for G R and Ga after multiplying with U(y, x) 
and taking the Wigner transform. Note that the right hand sides of eqs. (|2.46| ) define two 
new Wigner functions, p(k,X) and T(k,X), which are real quantities (cf. the discussion 
after eq. (|2.26| )) and can be seen as off-equilibrium generalizations of the corresponding 



spectral densities in equilibrium (recall eq. (|2.13|) ). In terms of these functions we have, 
for instance: 

G R (k,X) = r^k P^X) ; t R (k,X) = - r^ l ik>0 : k,X) ■ (2.47) 

j-oo 2ir k — Kq — ii] j-oo 2tt k — kq — irj 



Eq. ( |2.45| ) also involves: 

2ReG R = G R + G A , 2Re± R (k,X) = ± R + t A . (2.48) 



Further manipulations allow us to put eq. ( [2.45|) in the form: 



{2k-d k ReT)-{d x G < ) + (d x ReX) ■ (d k G<) 

(d k t<) ■ (d x ReG R ) + (d x t<) ■ (d k ReG R ) 

/ G > ± < - ± > G < ) , (2.49) 



where ReS = Ret R + S, and d x = d x - eF^(X)d^. It is interesting to note that 
the corresponding equation for a scalar field theory (like A0 4 ) can be obtained by simply 
replacing d^ by d x in the above equation 0. 

In equilibrium, both sides of eq. (|2.49|) are identically zero. This is obvious for the 



terms in the l.h.s., which involve the soft derivative d x or mean field insertions, and can 
be easily verified for the terms in the r.h.s. by using the KMS conditions for G eq and E eg 
(cf. eq. Q2.12j )). Thus, eq. ( |2.49|) is a transport equation which describes the space-time 



evolution of long- wavelength fluctuations in the average density of the charged particles. 
It holds to leading order in the gradient expansion (that is, up to terms involving at least 
two powers of the soft derivative), and to all orders in the interaction coupling strength. 

To conclude this section, note that, within the previous approximations (that is, 
up to terms involving at least two soft derivatives), the retarded propagator G R (k,X) 
satisfies an equation which is formally identical to that it obeys in equilibrium: 

(k 2 - t(X) - ± R (k, X)) G R (k, X) = -1 . (2.50) 



In order to obtain this equation, start with eq. ( |2.24| ) for G R (x, y) together with its 



conjugate equation where the differential operator acts on y; then, consider the sum of 
these two equations, and perform a gauge-invariant gradient expansion as above. In this 
expansion, all the terms involving one soft derivative d x cancel, and the same holds also 
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for the terms involving the soft mean field. From eqs. (|2.47 ) and ( |2.50| ), we deduce an 
expression for the off-equilibrium spectral density: 



r(k,x) 



p(k,X) = 2lmG R (k,X) = ■ 

(k 2 -ReE(k,X)) + [T(k,X)/2 

which will be useful in discussing the quasiparticle approximation below. 



2 ! 



(2.51) 



2.4 Mean-field and quasiparticle approximations 

In order to make progress with eq. ( |2.49| ) further approximations are needed. In particular, 
we shall use below perturbation theory to express the self-energies and S > in terms of 
the propagators G < and G > . As a first step, let us consider the mean field approximation 



in which the self energies E are neglected altogether. The equation (|2.49| ) reduces then 
to {k ■ dx)G K (k, X) = 0, or, more explicitly: 

k-d x -ek tl F^{X)d^)G < {k,X) = 0. (2.52) 

This equation describes the motion of independent particles in the mean field F^ v . In this 
approximation the spectral density remains the same as in the free theory in equilibrium, 
as obvious from eq. Q2.51|) : p(k,X) « po(k) = 2-ne{kQ)5{k 2 )N{ko). Accordingly, the 
solution to eq. ( |2.52[ ) can be written in the form: 

G < (k, X) = 2T,5(k 2 ){e(k )N + i}i, X) + 0(-&o)(l + N_(-k, X))}, (2.53) 



where the density matrices iV±(k, X) satisfy the Vlasov equation [p^j : 



(v-d x ± e(E + v x B) • V k )N±(k,X) = 0, (2.54) 

with f M = (l,v) and v = k is the velocity of the charged particle. The density matri- 
ces N±(k,X) may be given the interpretation of classical phase-space distributions for 
particles and antiparticles. In terms of them, the induced current is simply: 

j»{x) = ef-j^v*(N + QL,X)-N-QL,xj). (2.55) 

Going beyond the mean field approximation, we need to take into account the 
various effects of the self-energies. We shall concentrate here on a commonly used ap- 
proximation which consists in neglecting the broadening of the single-particle states when 
computing the collision terms, an approximation which we refer to as the "quasiparti- 
cle approximation". Indeed, the interaction rate T in eq. ( [2.5 Ij ) is of higher order in e 
(specifically, T ~ (e 2 ln(l/e))T 2 , as we shall see below), so we can use the mean field 



15 



spectral density, p(k,X) ~ Po(k), to estimate the collision terms. At the same time, we 
shall ignore the self-energy terms in the l.h.s. of eq. ([2.49 ). That this is consistent can be 



verified by power counting (we shall do this explicitly for the QCD case, in Sec. 3.5); it is 
also physically motivated from the fact that the role of these terms is to account for the 
difference between p(k,X) and po(k) in the transport equation for G < (fc,X) 0. Thus, in 
the quasiparticle approximation, the Wigner functions G < (k,X) and G > (k,X) preserve 
the same on-shell structure as in the mean field approximation, as displayed in eq. ( |2.53|) . 

We thus end up with the following kinetic equation: 

2{k-d x - ek ll F llv (X)dl) G < {k,X) = - - ± > G < ) , (2.56) 

where, in line with the weak coupling expansion, we choose the self energies S < and S > so 
as to reproduce the one-photon-exchange scattering in Fig. [I] (Born approximation). As 
we shall see in the next section, this generates a collision term of the standard Boltzmann 
form. Note also that eq. ( |2.56| ) in independent upon the choice of the path 7 in eq. ( |2.28| ); 
indeed, the terms which were explicitly path-dependent in eq. ( |2.49| ) have disappeared in 
the approximations leading to eq. ( (2.56j ). Moreover, we shall verify shortly that, to the 
order of interest, the self-energies S < and S > are path- independent as well. 

In computing transport coefficients like viscosities or electric conductivity (see, e.g., 



Refs. ||, [14], [T|, |34], f43|), or the quasiparticle damping rate [jy], g], it is only necessary to 
consider small off-equilibrium deviations, so that the linearized version of eq. ( |2.56| ) can 
be used. We then write, e.g., G< = Gf q + 5G< and £< = E< + <5£< (with 5G < G eq , 
<5E -C S eg ), and linearize the collision term with respect to the small fluctuations 5G and 
5±: 

C{k,X) = -(G > {k,X)± < {k,X) -£ > (k,X)G < (k,X)) 

£< 5G> - £> 5G<) + (5±> G< - 5±< G>). (2.57) 

In the quasiparticle approximation, we further have G> ~ Gq , Gf ~ Gq and 5G < ~ 
5G> = 5G (since G>(k,X) - G<(k,X) = p (k) = G$(k) - G$(k)). Then, the linearized 
collision term takes the form: 



C(k, X) ~ -T eg (k) 5G(k, X) + (<5£> G< - 5S< G>) , (2.58) 

where we have isolated the damping rate in equilibrium (cf. eq. (|2.46|) ): 

T eq {k) = S<(/c)-E>(A;). (2.59) 

Note that the quasiparticle approximation is not a self-consistent approximation, but it 
is in line with the weak coupling expansion: the collision term generates a width which is 
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k k-q 

Figure 3: The two-loop self-energy diagram which describes collisions in the Born approx- 
imation. The wavy lines denote free, equilibrium, photon propagators. The other lines 
are off-equilibrium scalar propagators. 

not included in the spectral densities which are used to estimate it; however, the neglected 
terms are of higher order than those we have kept. 

The most direct application of the formula above is the calculation of the quasi- 
particle damping rate PB] , |TT| , |3|. To this aim, we consider a specific off-equilibrium 
deviation which is obtained by adding, at t = 0, a particle with momentum p and energy 
Po = e p = |p| to a plasma initially in equilibrium. Since, for a large system, this is a 
small perturbation, we can neglect all mean field effects and assume N(p, t) to be only a 
function of time (here, N(p,t) = N + (p,t)\ cf. eq. ( 2.53Q ). Moreover, for momenta k ^ p, 



the distribution function does not change appreciably from the equilibrium value N(Ek), 
so that, to leading order in the external perturbation, we can ignore the off-equilibrium 
fluctuations of the self-energies: 5±< w 5±> « 0. Then, eqs. ffT53D , ( p36|) and ( ggg ) 
yield a very simple equation for the fluctuation 5N(p,t) (with T(p) = T eq (p = s p ,p)) : 

2e p pN(p,t) = -T(p)5N(p,t), (2.60) 

whose solution shows exponential attenuation in time: 

8N(p,t) = 5N{p, 0) e" 27(p)t . (2.61) 

The quasiparticle damping rate 7 is here conventionally defined as ^(p) = r(p)/4e p P]. 
(This simple picture is actually complicated by infrared effects to be discussed in Sec. 2.6 
© I§.) 
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Figure 4: The effective one-loop self-energy which describes collisions in the resummed 
Born approximation. The blob on the photon line denotes the resummation of the one- 
loop polarization tensor (cf. Fig. @). 



Figure 5: One of the multi-loop diagrams (containing n bubble insertions along the photon 
line) which is included in the effective one-loop diagram in Fig. |j. 

2.5 The collision terms 

We now turn to the calculation of the collisional self-energy corresponding to Fig. p]. As 
it is well known (and will be verified later), the corresponding transport cross section 
is dominated by relatively hard momentum transfers, eT ^ q ^ T. When the photon 
momentum is hard, q ~ T, the process in Fig. [I] is described by the two-loop self-energy 
depicted in Fig. |3] in which all the lines are hard, and the scalar propagators are to be 
understood as off-equilibrium propagators. The photon propagators, on the other hand, 
are just free propagators in equilibrium. At soft momenta q ~ eT, the relevant self- 
energy is given by the effective one-loop diagram in Fig. |] in which both the internal lines 
denote off-equilibrium propagators; the scalar line is hard, while the photon line is soft 
and dressed by the off-equilibrium polarization tensor in the one-loop approximation (this 
is denoted by a blob). That is, the diagram in Fig. [| involves an infinite series of bubble 
insertions along the photon line, as illustrated in Fig. |5|. 

It is furthermore convenient to recognize that the diagram in Fig. || is one of the 
family of diagrams displayed in Fig. |5]; thus, we can use the effective one-loop self-energy 




k 



k-q 
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in Fig. |] to describe the collision in Fig |T| for all photon momenta, which we shall do in 
what follows. To evaluate this diagram, we need the vertex coupling the photon to the 
scalar field in the presence of the classical background field A^. This can be read off the 
Lagrangian: 



C = {pi - iey) f (D" + iea M ) 



(D^t-tfiD^) +eV^V, (2-62) 



2 n n»Aj, 



where _D M = <9 M + ieA^. There are two relevant vertices: —2ieft a^D^cj) and —ie 
The self-energy reads then: 



'i-i- 



E(x,y) = e 2 ^{D^G(x,y))v, u (x,y) + 2{D^y))(d;V flu (x,y)] 

+ 2(D^G(x,y))(dZV, u (x,y))+G(x,y)(d;d»V^(x,y))\. (2.63) 



Here, V^ix.y) is the off-equilibrium photon propagator, to be constructed shortly. By 
appropriately choosing x and y along the contour, we get expressions for both and 
S < . For instance, S > (x, y) will involve G > (x,y) and D > (x,y), etc. Below, to simplify the 
notations, the upper indices > and < will be often omitted. 

We need then to evaluate the gauge invariant self energy T,(x,y) = U(y,x)Ti(x,y). 
In doing that, we meet terms like: 

U(y, x) D% G(x, y) = U(y, x) D% (u(x, y)G(x, y)) . (2.64) 

Performing the gradient expansion of such a term, one gets: 

U(y,x)D% (u{x,y)G{x,y) 
« e ies - A{x) (d% + ieA"(X)) {e~ ies - A{x) G(s,X)} 
« d»G(s,X). (2.65) 

Note that, in the above manipulations, we have used the simple approximation ( |2.31| ) for 
U(x,y), which makes the final result independent of the choice the path in the Wilson 
line. The same holds for all the other results in this section. 

Similarly, we get dyT> fll/ (x,y) ~ — d u s V ^ u {s , X) , so that the expression 
D%G(x,y)dyV flu (x,y) becomes simply —d£G(s, X)dg / D fll/ (s, X). Proceeding in the same 
way for the other terms, and performing the Wigner transform, we get: 

±(k, X) « e 2 J j0p ( 2 p - tf) (2F - <f ) G(k - q, X)V„ v (q, X). (2.66) 

This is formally the same expression as in equilibrium, except for the fact that G(k — q, X) 
and V^ v (q, X) are off-equilibrium propagators and fc M has to be interpreted as the kinetic 
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(a) (b) 

Figure 6: One-loop contributions to the photon self-energy in SQED. All the internal lines 
are off-equilibrium propagators. 

momentum. (Note that the photon propagator V^ v does not need a special treatment 
since it is invariant under the gauge transformations of the background field.) 

The photon propagator D^ix, y) obeys Dyson- Schwinger equations similar to eqs. 
( ggp . Specifically, 

(gi»& - - U^) x V<{x,y) = J d 4 z (U R V< + U<V A Y p (x,y), (2.67) 
and similarly: 

(g»ud 2 - d,d u - fl„ v ) x V%(x,y) - J d"zU^(x,y) = 8>SP{x - y) , (2.68) 

where, to the order of interest, U^x, y) is given by the one-loop diagrams in Fig. |. 
That is, H^ u (x) is the tadpole contribution in Fig. ^.b, while the non-local self-energies 
Hr(x, y) and U < (x, y) are determined by the graph in Fig. |6].a. From the equations above, 
we deduce the following relation between T> K and Il < : 

V%(x,y) = - J d A Zl d 4 z 2 (v R (x,z 1 )Il < (z 1 ,z 2 )V A (z 2 ,y))^ (2.69) 

which becomes, after a gradient expansion, 

P< {q, X) » - (p R {q, X) U<(q, X) V A {q, X)) ^ (2.70) 

up to corrections of 0(dx/q)- Since, as we shall see shortly, the collision terms are 
saturated by momenta q ^ eT, the corrections to eq. ( |2.70| ) are of higher order in e 
provided dx ~ e 2 T . A similar relation holds between T> > (q, X) and n > (g, X). 

It should be observed here that a new scale is entering the gradient expansion. In 
most situations before, the soft derivative dx appeared in combinations such as s-dx with 
the magnitude of the non- locality fixed by thermal fluctuations: s ~ 1/T. In eq. (|2.69| ). 
however, the non-localities x — Z\ or z 2 — y are of order 1/q and may be interpreted as 
the range of the effective interaction between the colliding particles. Thus, the validity of 
the gradient expansion in this case relies on the range of this effective interaction being 
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2.22)- 



small compared to the scale of the inhomogeneities, as measured by c^ 1 . Now, the range 
of the effective interaction depends on the specific transport processes one is looking at. 
In most cases, and as a result of cancellations to be exhibited in the next subsection, this 
range is typically of order 1/T (and marginally 1/eT) so that the gradient expansion is 
indeed valid to calculate transport coefficients to leading order in e already for processes 
taking place on a scale l/e 2 T. 

To construct the photon self-energy out of equilibrium, we use the interaction ver- 
tices from eq. (|2.62 ) and obtain, to the order of interest, 

IV(X) = -2g, u e 2 J *L6<(p,X), (2.71) 

and (compare to eq. ( |2.66| )) : 

U> u (q,X) « e 2 J -j^- (2p, + q f ,){2p v + q u )6>{p + q,X)G<{p,X), (2.72) 

together with a similar expression for H < (q, X) which involves G < (p+q, X) and G > {p, X). 
These expressions are gauge invariant, as expected. The tadpole piece ( |2.71|) enters the 
calculation of the retarded propagator T>n(q, X), which is related to the self-energy fl(X) + 
Uii(q, X) by the same equation as in equilibrium (cf. eq. ( |2.50| )). 

By collecting the previous results, we finally obtain the following collision term: 



d 4 p r d A q 
{2nyl (2tt) 



C(k,X) = - \M pk ^\ 



x {G<{k, X)G<(p, X)G>{k', X)G > {p', X) - G > (k, X)G > {p, X)G < {k', X)G < {p', X)} , 

(2.73) 

where p' = p + q, k! = k — q, and Ai p k-> p >k' is the scattering matrix element: 

\M pk ^p'k'\ 2 = e\k + k'Yip + p'Y{k + k') a (p + p'f[D R (q, X)]^[V A (q, X)U ■ (2.74) 

This collision term has the standard Boltzmann structure, with a gain term and a loss 
term. To be in line with the previous approximations, this must be evaluated with the 
Wigner functions G < and G > in the quasiparticle approximation, i.e., G < (k,X) has the 
on-shell structure exhibited in eq. (|2.53|) , while G > (fc, X) reads similarly : 

G > (k,X) = 2n5(k 2 ){6(k )(l + N + (k, X)) + 6(-k )N_(-k, X)}. (2.75) 

Thus, the four energy variables k , p , k' and p' in eq. ( p.73j ) are always on shell (e.g., 
|fc | = Ek = |k|), but they can be either positive (k = e fc ) or negative (k = — corre- 
sponding to particles and antiparticles, respectively. Choose k > for defmitness; then 



21 



p p -p -p 

► 7 ► < 7 < 



(a) (b) (c) 

Figure 7: Scattering processes described by the collision term in eq. (|2.73|) : (a) particle- 
particle scattering; (b) particle-antiparticle scattering (t-channel); (c) particle-antiparticle 
annihilation (s-channel). In Fig. (c), the virtual photon is always hard, so it needs no 
resummation. 



( |2.73|) describes the t-channel particle-particle scattering depicted as Fig. [I] (when all the 



energy variables are positive), but also the t-channel particle-antiparticle scattering (when 
k' is positive, but po and p' Q are both negative), and the particle-antiparticle anihilation 
(or s-channel scattering: p and k' negative, and p' positive). These various processes 
are illustrated in Fig. [7|. 



2.6 Quasiparticle lifetimes vs. relaxation times 

At the end of Sec. 2.4, we have seen that the collision term fl2.58| ) yields the quasiparticle 
lifetime r ~ I/7, which is dominated by soft momentum transfers, q ^ eT, and is 
typically r ~ l/(e 2 Tln(l/e)) |3l], |32| . For transport phenomena, however, it is well 
known 13, |HJ that the Abelian collision term ( |2.73| ) is saturated by relatively large 
momentum transfers eT ^ q ^ T, and the typical relaxation time for off-equilibrium 
perturbations is r tr ~ l/(e 4 Tln(l/e)). The fact that r tr ^> r is due to specific infrared 
cancellations in the collision term, that we shall discuss now. Note that in QCD the colour 
algebra prohibits similar cancellations in the calculation of colour relaxation processes, as 
we shall see in Sees. 3.8 and 4 below. 

Consider then the linearized version of the collision term, as given by eq. fl2.58| ). It 
is convenient to define (with (3 = 1/T) : 

diV 

5G(k,X) = Po (k)5N(k,X) = -p (k)W(k,X) — 

clfc 

= P Po (k)W(k,X)N(k )[l + N(k )} (2.76) 

where (cf. eq. Q2.75p ) : 

5N(k,X) = 9(k )5N+(k,X)) - 9(-k )5N„(-k,X). (2.77) 
The function W(k,X) describes the local distorsion of the momentum distribution, as 
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may be seen from the following equation: 

dN 

N(k,X) = N(ko) - W(k,X) — « N(ko-W(k,X)). (2.78) 

dfco 

In terms of these new functions, the linearized collision term takes a particularly simple 
form, 

C(k,X) « -Pp (k)JdT\M\ 2 N(k )N(p )[l + N(k' )][l + N(p' )] 

x [W{k, X) + W(p, X) - W(k', X) - W(p', X)}, (2.79) 
with the following notation for the phase-space integral: 

/ dT = / / Po(P)Po(P + - q). (2.80) 

The matrix element \Ai\ 2 in eq. ( [2.79| ) is to be computed with the equilibrium retarded 
and advanced photon propagators V^A^q) (cf. eq. ( |2.74|) ). 

Following eq. ( |2.58| ), we identify the damping rate T(k) = T eg (k) as the coefficient of 
5G(k, X) in the r.h.s. of eq. ( |2.79| ), that is, as the term involving the fluctuation W(k, X) : 

C x {k,X) = -Ppo(k)W(k,X)Jdr\M\ 2 N(k )N(p )[l + N(k' )][l + N(p' )] 

= -T(k)8G(k,X). (2.81) 

One can verify that within the present approximation the above expression of T satisfies 
indeed eq. ( |2.59| ). 

Moreover, in eq. ( p.81|) , T(k) must be evaluated on the tree-level mass-shell (i.e., at 
l^o I = £fc), since it is multiplied by the on-shell fluctuation 8G(k, X). This determines the 
quasiparticle damping rate, 7 = (l/4k)T(ko = k), which is, however, well known to be 
infrared divergent in the present approximation |28], |30], [31], [32| . Specifically, the leading 
contribution to 7 comes from soft momenta exchange q ^ eT in the t-channel collisions 
in Figs, ^.a and b. To evaluate this contribution, we can neglect go next to po and ko in 
the thermal distributions in eq. (|2.81|) , and get: 

T(k)~ JdT\M\ 2 N(p )[l + N(p )}. (2.82) 

To the order of interest, we need the resummed photon propagator in the "hard thermal 
loop" approximation 0,0, to be denoted as \q). This yields (for q k, p): 

\M\ 2 ~ 16e 4 44 \* V M) + (q x v) ■ (q x v') *V t {q)\ 2 . (2.83) 

where *T>\ and *T> t are the longitudinal (or electric) and the transverse (or magnetic) 
components of the retarded propagator, with the following IR behaviour (below, mo is 
the Debye mass, m 2 D = e 2 T 2 /3) 0, g : 

*P,(g Q -> 0, q) ~ ~ X 2 , *P t (g Q « g) ~ ^ . 1 - 2 . (2.84) 

g^ + m|) q z — t {Ttqo/Aq) 



23 



Because of Debye screening, the electric contribution to the damping rate 7/ is finite and 
of order e 4 T 3 /m 2 D = 0(e 2 T). In the magnetic sector, the dynamical (go 7^ 0) screening 
JT|, H is not enough to make finite 74, which remains logarithmically divergent (see Ref. 



32] for more details): 



6 4 T 3 - .g ^ 1 

It - ~ttt~ I dg / — 



24 J J- q 2tt g 4 + (vrm^go^g) 2 
e 2 T /" m o dg e 2 T m D 



, , , In—, (2.85) 

An q An ^ 

where /1 is an IR cutoff and we have retained only the dominant, logarithmically divergent, 

contribution. The remaining IR divergence in eq. ( 2.85|) is associated to the unscreened 



static magnetic interactions. Since the latter have an infinite range, one may worry that 
the gradient expansion may become invalid in the calculation of the damping rate (cf. the 
remark after eq. ( |2.7(J| )). Recall, however, that in the calculation of 7, the particles with 
which the quasiparticle interacts are in equilibrium, and constitute therefore a uniform 
background (cf. the discussion at the end of Sec. 2.4.) Thus, the question of the relative 
sizes of the range of the interaction and that of the space-time inhomogeneities is not 
an issue here. Rather, the IR divergence in eq. Q2.85 ) is an artifact of the perturbative 



expansion and can be eliminated by a specific resummation [^] which goes beyond the 
approximations performed in deriving the Boltzmann equation (see however Ref. 



The IR problem of the damping rate does not show up in the calculation of the 
transport coefficients, because the IR contribution to T, the first term in the r.h.s. of 
eq. ( |2.79| ), is actually compensated by a similar contribution to the third term, involving 
W(k',X) : indeed, for soft q, W(k',X) = W(k - q,X) « W(k,X), so that the first 
and third terms in eq. (|2.79|) cancel each other. As we shall see in Sec. 4, this can 



be understood as a cancellation between self-energy and vertex corrections in ordinary 
Feynman graphs. A similar cancellation occurs between the other two terms in eq. ( |2.79|) , 
namely W(p, X) and W(p f , X). Thus, in order to see the leading IR (g <C T) behaviour 
of the full integrand in eq. ( |2.79| ), one has to expand W(k',X) and W(p',X) to higher 
orders in q. This generates extra factors of q which remove the most severe IR divergences 
in the collision integral. As a result, the typical rate involved in the calculation of the 
transport coefficients is T tr ~ e 4 Tln(l/e), where the logarithm originates from screening 
effects at the scale eT. 

Of course, the simple arguments above are only good enough to provide an order- 
of-magnitude estimate for the transport relaxation times. In order to compute transport 
coefficients, one has to solve the Boltzmann equation ( ]2.56|) with the linearized collision 



term (|2.79|) , which is generally complicated. Explicit solutions can be found, e.g., by 
using specific Ansatze for the unknown function W(k,X), or by variational methods. 
Some calculations of this kind can be found in Refs. 0, |14], ITBT |33l |34|. 
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3 Boltzmann equation for hot QCD 



We now come to the case of the high temperature Yang-Mills plasma. As mentioned in the 
Introduction, we are interested in the regime of ultrasoft colour excitations propagating 
on a typical scale A ~ l/g 2 T. (More precisely, the spatial gradients of the fields are of 
order g 2 T, but their time derivatives can be even softer, i.e., of order g A T.) The relevant 
response function is the induced colour current, which we shall eventually express in the 
form (with = (1,1s./ k)): 

3 ;(X) = ^J-0- 3 v,Tr (T«6N(k,X)). (3.1) 

where 5N ab (k, X) is a density matrix in colour space. The overall factor 2 stands for the 
two transverse polarizations. 

The density matrix 5N ab (k, X) is a functional of the average fields and must 
transform covariantly under the gauge transformations of the latter. That is, under the 
gauge transformation (A^ = A a ^T a ): 

A"(X) — ► h(X) (a^X) - ~ d^J rf(X), (3.2) 

where h(x) = exp(i8 a (x)T a ), we must have: 

5N ab (k,X) — h a - a (X)5N ai (k,X)hl(X). (3.3) 

Indeed, this ensures that j^(X) transforms as a colour vector: j£ — > habjtf, or, in matrix 
notations, 

f(X) = j%T a — h(X)f(X)h\X). (3.4) 

(This should be contrasted with the Abelian case, where both the current and the distri- 
bution function are gauge invariant.) 

The covariance of the density matrix SN a b()s, X) should result from a corresponding 
property of the off-equilibrium gluon propagator G^(x,y) = (Ta%(x)a%(y)) from which 
it originates. However this propagator depends not only upon the choice of a gauge 
for the average field A%, but also on the gauge- fixing condition for the fluctuating field 
a%. With a generic gauge fixing, G^(x,y) transforms in a complicated way under the 
gauge transformations of A%. The situation becomes simpler when one uses the so-called 



"background field gauge" to be introduced in the next subsection |45], |46|. Then the gauge 
fixing term is covariant under the gauge transformations of the average field A%, and the 
gluon propagator G(x, y) can be turned into a covariant quantity by attaching Wilson 
lines in x and y. We shall then be able to maintain explicit gauge symmetry with respect 
to the background field at each step of our calculation. 
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3.1 The background field gauge 



In this method, one splits the gauge field into a classical background field A a ^ to be later 
identified with the average field, and a fluctuating quantum field a". The generating 
functional of Green's functions is written as: 

Z\j; A}= J VaV(V( e iSpp[a ^' A]+i Ic d4x ^< , (3.5) 

with the Fadeev-Popov action: 

S FP [^ C, C; A] = J c d 4 x{- \{F%,[A + a]) 2 + ±.(D t [A}af + ( a (p i [A\D i [A + a})j b }, 



(3.6) 



where D^[A + a] = <9 M + ig(A^ + a M ) is the covariant derivative for the total field A^ + a M , 
and F® u [A + a] is the respective field strength tensor. Furthermore, the gauge-fixing term 
{l/2\)(Di[A]a 1 ) 2 , which is of the Coulomb type, is manifestly covariant with respect to the 
gauge transformations of the background gauge field A^. Accordingly, the exponential in 
eq. ( |3.6|) is invariant with respect to the following transformations (with matrix notations: 
h(x) = exp(i6 a (x)T a ), a M = a^T b , ( = ( a T a , etc.): 

Apu -> h{Ap - (i/g)d fl )h\ j M -> hj^h), 

a„ -> Zia^, C -> K^, C -> ^ f Cfr- (3-7) 

(Note the homogeneous transformations of the quantum gauge fields (a M ) and ghost fields 
(C, C) i n ^e equations above.) Because of this symmetry, the generating functional Z[j; A] 
is invariant under the normal gauge transformations of its arguments, given by the first 



line of eq. (|3.7| ). Then, the gluon Green's functions, derived from Z[j; A] by differentiation 
with respect to j°, are gauge covariant under the same transformations. 

The physical Green's functions are obtained by identifying the total average field 
to the background field. This implies: 

, i, .. SlnZlj; A] 

«M> - = 0.8) 

which determines a functional relation between the external current and the average field; 
we write this as j = j[A]. Then, the 2-point function is obtained as: 



Gj(x,») = (T c a*(xK(y)> - ' )JluZLy: l! 



(3.9) 



Under the gauge transformations ( |3.7| ) of A^, it transforms covariantly: 

G%(x,y) - h a - a (x)G%(x,y)hl(y). (3.10) 
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The ghost propagator, 



A ab (x,y) = (T c ( a (x)( b (y)), (3.11) 

has the same transformation property. Similar covariance properties hold for the higher 
point Green's functions, and for the various self-energies. Note that, in practice, we shall 



never have to solve the implicit eq. ( |3.8|) for j[A], since we shall be able to impose the 
condition (a^x)) = directly on the equations of motion for the Green's functions. 

In deriving the Boltzmann equation satisfied by 5N(k, X), it will be convenient to 
use the Coulomb gauge, which offers the most direct description of the physical degrees of 
freedom: in this gauge, the (hard) propagating modes are entirely contained in the trans- 
verse components of the spatial gluon propagator Gij(x,y), so that the density matrix 
5iV(k, X) is simply the gauge- covariant Wigner transform of Gij(x,y) (see below). (In 
other gauges — like the "covariant" ones with gauge-fixing term (l/2\)(D^[A]a f1 ) 2 — the 
physical, transverse degrees of freedom are mixed in all the components of the gluon prop- 
agator G^y. In this case, the density matrix 5N(h, X) involves a linear combination of the 
Wigner functions of the gluons and the ghosts, and it is only this particular combination 
which is gauge- fixing independent [[|. The intermediate calculations are cumbersome, 
and the explicit proof of the gauge-fixing independence is quite non-trivial already at the 
mean field approximation — or "hard thermal loop" — level |2£|, f|].) 

In what follows we shall mostly use the strict Coulomb gauge condition, namely: 

A[^K = 0- ( 3 - 12 ) 
In this gauge, all the non-equilibrium Green's functions are transverse, that is: 

Di[A]G iu (x,y) = 0, (3.13) 

and similarly for the higher point functions. The only non-trivial components of the free 
retarded gluon propagator are: 

G®(k) = -± 4 0) W = -^^- (3-14) 

That is, the electric gluon is static, and the same is also true for the Coulomb ghost: 
A^(k) = 1/k 2 . Accordingly (with G$(k) and Gft(Jfe) as defined in eq. (HI)), 

Gf^\k) = {5ij - hk) Go(k), Gl i0 \k) = {5ij - hk) G£{k), (3.15) 

while all the other components are zero. 
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3.2 Equations of motion 

The equations of motion for the average field A p read: 

(D»F^) a (x) = j;(x). (3.16) 

Here and in what follows, or denote the covariant derivative or the field strength 
tensor associated to the background field A ^. The induced colour current j°l{x) involves 
the off-equilibrium 2-point functions^ for gluons and ghosts: 



f a {x) = ig^T a {T^DlG< u {x,y) + A<(x,y)(D^} 



(3.17) 

y=x 



We have used here the notation: 

r ^pA = 2g^g px - g pp g vX - g px g vp . (3.18) 

Furthermore, -D^[^4] = d — igA a T a , and the derivative d acts on the function on its left. 
The Kadanoff-Baym equations for the gluon 2-point functions read (cf. Sec. 2.1) : 

(gPD 2 -D^ + 2igF;) x G< u (x,y) = 

Jd 4 .{ V E^x,z)G<(z )2 /) + Z< p {x,z)G p A x {z,y)g x „}, (3.19) 

and 



G< p (x,y)(g pv (rt) -DlDl + 2igF pL 

Jd 4 z{g, x G x R p (x,z)E< u (z,y) + G< p {x,z)E pX {z,y)g Xu }, (3.20) 



together with the gauge fixing conditions (cf. eq. (|3.13j )): 

DiG iu (x,y) = 0, G^(x,y)D^ = 0. (3.21) 
In deriving these equations, we have used symmetry properties like: 

G>*(x,y) = G< b ;(y,x), G R ab ^(x,y) = Gfafax), (3.22) 

and similarly for the self-energies. 

In the following developments, we shall often omit the upperscripts > and < on the 
2-point functions, and indicate them only when necessary, e.g., on the final equations. 



e There is also a contribution to the current from the gluon 3-point function which, however, starts at 
two-hard- loop level and is thus negligible for what follows || . 
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3.3 Gauge-covariant Wigner functions 



Let G a b(x, y) denote any of the 2-point functions, and G a b(k, X) the corresponding Wigner 
function, defined as in eq. Q2.26Q . Unlike G a b(x,y), which is separately gauge-covariant 
at x and y (cf. eq. ( 3.10Q ), its Wigner transform G a b(k, X) is not covariant. However, 



following what we did for SQED, we can construct the following function (cf. eq. ( j2.29| )): 

G ab (s,X) = U a - a (x,X+^)G- ab (x + ^X-^)U bb (X-^,x), (3.23) 

where U(x, y) is the non-Abelian parallel transporter, also referred to as a Wilson line 
{A, = A^T a ) : 

U{x, y) = P exp j-ip J dz'A^z)} . (3.24) 

As in the Abelian case, the path 7 is arbitrary (see the discussion before eq. ( |2.31| )). 
Under the gauge transformations of A^, the Wilson line ( (3.241) transforms as (in matrix 
notations): 

U(x,y)^h(x)U(x,y)h^y), (3.25) 

so that the function ( |3.23| ) is indeed gauge-covariant at X for any given s: 

G{s,X) — > h(X)G(s,X)h i (X). (3.26) 

Correspondingly, its Wigner transform G a b(k,X) transforms covariantly as well: For any 
given k, G(k, X) — ► h(X) G(k, X) h\X). 

In principle, the equations of motion for G(s, X) follow from the equations of motion 
(l£igD-(l£2TD f o r G{x,y) by replacing G{x,y) by (cf. eq. (ggg)): 



G(x, y) = U{x, X) G(s, X) U{X, y). (3.27) 

However, in contrast to what we did for SQED, in the non Abelian case we have to proceed 
to a linearisation in order to preserve the consistency of the expansion in powers of g. 
Recall indeed that the mean fields are supposed to be weak and slowly varying, such 
that dx ~ gA <C T. (The ultrasoft covariant derivative is of the order Dx = 0(g 2 T), 
but the simplifications we are refering to hold already when Dx = O(gT) [§, |], HO 
For such soft background fields the function G(s,X) remains strongly peaked at s = 0, 
and vanishes when s ^ 1/T. Over such a short scale, the mean field A^ does not vary 
significantly. Furthermore, for s ^ gs ■ A 1 since gA^ T. We can then expand 
the Wilson lines in eq. ( 3.23|) in powers of g and get, to leading non-trivial order: 



U ab {x,y) ~ 5 ab - ig(s- A ab (X)). (3.28) 
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This should be compared to eq. ( [2.3 1| ) in SQED: both expressions hold to leading order 



in an expansion in soft gradients, but in the non-Abelian expression ( 3.28 ) we have also 
performed an expansion in powers of the gauge field. In what follows, we will never need 
to go beyond the simple approximation (|3.28| ). 

Similarly, we shall see that the off-equilibrium fluctuations SG = G — G eq are per- 
turbatively small: SG ~ (D x /T)G eq ~ g 2 G eq . Thus, by writing: 

G = G eq + 5G, G = G eq + 8G, (3.29) 

in eq. ( |3.23| ), and recalling that G ab = 5 ab G eq , we can easily obtain the following relation 
between 5G and SG, valid to leading order in g: 

66(8, X) ~ 6G(x,y) + ig(s-A(X))G eq (s), (3.30) 

or, equivalently: 

6G(k, X) ~ 6G(k, X) + g(A(X) ■ d k )G eq (k). (3.31) 

Note that both terms in the r.h.s. of eq. ( |3.30| ) or (|3.31|) are of the same order, namely of 
0((.D x /T)G eg) . On the other hand, the terms which have been neglected in going from 
eq. ( |3.23| ) to eq. ( |3.30| ) are down by, at least, one more power of D x /T. 

Consider now a term like D£G(x,y) which appears in eqs. (|3.19|) -( |3.21|) . Clearly, 



such a term transforms in the same way as G(x,y), so it can be treated in a similar way 
(cf. eq. ( |3.23|) ). Then, we can write: 



6(D^G(x,y))=D^G(x,y)-d^G eq (s) ~ 8*86(3, X) - ig(s ■ A(X))d»G eq (s), (3.32) 



which parallels eq. ( |3.30|) . In particular, since the equilibrium gluon Wigner function 



is transverse, diG™ = 0, eqs. ( |3.21|) and ( |3.32|) show that the gauge- covariant Wigner 
function is transverse as well: 

8 i s 6 iv (s,X) = 0, or VGU^X) = 0. (3.33) 



Finally, we have to express the induced current (|3.17 ) in terms of the gauge-covariant 



Wigner functions. Since it vanishes in equilibrium, it involves only the off-equilibrium 
deviations of the Wigner functions of the gluons and the ghosts. We have: 

j;(X) = g J 0^TrT a {-k,5Gr(k,X) + 8G< u (k,X)k»-k fl 8A<(k,X)}, (3.34) 



where the following property has been used (cf. eq. ( |3.32| )): 



y=x 



D%G(x,y) =d<?G(s,X) . (3.35) 



s=0 
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Like ( |3.17| ), eq. ( |3.34| ) holds in an arbitrary gauge. In Coulomb's gauge it can be further 



simplified: as we shall see in the next section, only the transverse fluctuations 5Gij(k, X) 
(5ij — kikj)5G(k, X) matter for the calculation of j 74 , so that: 

d 4 k 



f a {X) = 2gJ — k »Tr{T a 6G<(k,X)}. (3.36) 



3.4 The non-Abelian Vlasov equation 



In this section, we shall study eqs. ( |3.19j )-( |3.21| ) in the limit where the all the terms 



involving self-energies can be neglected. As in the case of SQED, this amounts to a mean 
field approximation in which the hard gluons are allowed to scatter on the average colour 
fields A^, but not among themselves. The resulting equations are: 

(g^D 2 - D^D" + 2igF f f) x G pu (x,y) = 0, 
G^y)(g pv {D^f-DlDl + 2igF f)V ) y = 0, (3.37) 

where G denotes either one of the functions G > or G < . The outcome of the present 
subsection is the Vlasov equation for the gluon density matrix. The derivation is not new 
0], except for the use of the Coulomb gauge. However, since this involves manipulations 
which will be essential for the evaluation of the collision terms, we present it in detail. 

The equations (|3.37|) involve hidden powers of g, associated with the soft inho- 
mogeneities (dx ~ 9 2 T) and with the amplitudes of the mean fields (A ~ gT and 
gF^v ~ g 4 T 2 ). The purpose of the covariant gradient expansion is precisely to isolate 
all the terms of leading order in g. (Actually, all the manipulations in this subsection 
apply already for inhomogeneities at the scale gT, when dx ~ gA ~ gT and gF pu ~ g 2 T 2 



As in Sec. 2.3, we start by considering the difference of the two equations ( |3.37|) . 



Let us look at the first term in the l.h.s. of this difference equation, which we denote as: 

E(x,y) = D 2 x G(x,y) - G(x, ^(Dj) 2 , (3.38) 

where D 2 and (-D^) 2 are given by eq. (|2.35|) , except that the derivatives in (D^ y ) 2 are now 
understood to act on their left. (Minkowski indices are omitted to simplify the notations; 
they will be reestablished when needed.) Proceeding as in Sec. 2.3, and paying attention 
to the colour algebra, we obtain: 



E(s,X) = 2d s -d x G + 2ig[A fl (X),d»G]+ig{A tl (X),d x G}+ig{(s-d x )A tl ,d2G} 
+ig {(d x ■ A),G} - g 2 [A 2 (X),G] - ^ {(s ■ d x )A 2 ,G} + (3.39) 
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where the right parantheses (the braces) denote commutators (anticommutators) of colour 
matrices, and the dots stand for terms which involve at least two soft derivatives d x . 

At this point, we use the fact that A ~ gT and 6G = G — G eq ~ g 2 G eq (as will be 
verified a posteriori), with G eq ~ G^ in the mean field approximation. To leading order 
in g, eq. (|3.39|) then simplifies to: 

E(s, X) « 2(d s ■ d x )6G + 2ig [A„ d^JG] + 2ig{s ■ d x )A fl (d»G®) + 2ig{d x ■ A)G<® , 

(3.40) 

where all the terms are of order g 4 T 2 Go. Taking now the Wigner transform, we get: 

E(k,X) « 2[k-D x , 5G(k,X)]+2gk»(d x A^X))d u G {0 \k), (3.41) 

where G(k, X) is the ordinary Wigner transform of G(x, y), defined as in eq. fl2.26p . This 
can be rewritten in a gauge-covariant form by replacing 6G = 6G — g(A ■ dk)G^ (cf. 
eq. (TO ): 

E^(k,X) « 2[k-D x ,6G^(k,X)] -2gk a F af3 (X)d^G^(k), (3.42) 

where the Minkowski indiced have been reintroduced. 

We return now to eqs. ( |3.37| ). Since we are mainly interested in the transverse gluon 
Wigner function SGtj(k, X), let us focus on the components \i = i and v = j: 



D\G l3 - DfD* G oj + 2tgF/(x)G 



pj 



G tJ (Dl) -G lQ Dl y D] y + 2igG ip F»M = 0. 



(3.43) 



(In writing these equations, we have also used the gauge-fixing constraint (|3.21|) to simplify 
some terms.) When taking the difference of these equations, we first meet (cf. eq. Q3.42] )): 

, 2 



D x Gij — Ga ( D 



k ■ D x , 5GJ - 2gk a F a0 (X) d^G§\k). 



(3.44) 



Note the following identity, which will be useful later: 

k a F a/3 d /3 G\f(k) = k a F a ^[{5 l3 - kkj)G (k)] 

= fa - kk^F^Go - k a F Q i ^ii + W*-*Mih Gq . 



(3.45) 



The terms involving G;o and Gqj vanish in equilibrium, and remain small out of equi- 
librium, but nevertheless their conribution to eqs. ( |3.43| ) is non-negligible: Indeed, we 
shall verify shortly that Gio ~ (D x /T 2 )Gq, which is one order higher than the transverse 
fluctuations 5Gij ~ {D x /T)Gq. However, the hard derivatives multiplying do and Goj 
in eqs. (|3.43j) do not cancel in the difference of the two equations, in contrast to what 
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happens with the spatial components SGij. Specifically, <9 2 — <9 2 = 2d s ■ dx ~ TDx, while 



df&Q ~ d?d^ ~ T 2 . Therefore, the difference: 

D?D%G 0j - G i0 Dl y D} y « ^(fijfGoi - S^w) ~ (D X ) 2 G Q , (3.46) 

is of the same order as, e.g., (d s ■ dx)Gij. We thus have to evaluate these terms properly, 
which we shall do later, with the following results: 

G oj (k,X) » 2igF 0l 5lj ~J lkj G (k), G$(k,X) = G%(-k,X). (3.47) 

(The second equality above follows from the symmetry property ( |3.22j ).) The correspond- 
ing contribution to the kinetic equation for 8G{j(k,X) reads then: 

- (DfDSGy - G i0 Dl y D) y ) — -2gk°F ol ^ + 2k * k A G ^ k y (3 . 48) 

Finally, in the last terms in eqs. (|3.43|) — the terms involving the field strength tensor - 
we can replace F i p {x) ~ Ff{X) and G P j(k,X) « Gpj(k) = 8 p i(8ij — kikj)Go, to get: 

- 2ig(F u (X)G\f(s) - G^&F^X)) — > -2gG (k) (Fukkj + F/hh). (3.49) 

By using the identity ( |3.45| ), it is easy to recognize the role of the two contributions in 
eqs. (|3.48| ) and (|3.49|) : this is to cancel the non-transverse piece in the r.h.s. of eq. ( 3.44|) . 
Finally, 5Gij(k,X) satisfies the following kinetic equation: 

k ■ D x , SGyfcX)] - - k t k J )9k a F a( ,(X)d (3 G (k) = 0. (3.50) 

It is transverse, as anticipated: 

8Gij(k,X) = (Sij - kikj)8G(k,X), (3.51) 

with the new function 5G(k,X) satisfying: 

k ■ D x , 8G(k,X)] = gk a F^(X)d (3 G (k). (3.52) 

Since k ~ T, D x ~ g 2 T and gF afi ~ (D x ) 2 ~ # 4 T 2 , eq. ( g5g ) implies ~ (D X /T)G ~ 
g 2 Go, as anticipated. 

Eq. (|3.50|) is the main result of this subsection. In order to complete its proof, we 
still have to justify eq. fl3.47| ) for Goj. To this aim, we shall consider the first eq. ( |3.37| ) 
with \i = and v — j. This reads 

D 2 G 0J + 2igF i(x)Gij(x,y) = 0. (3.53) 

To the order of interest, D 2 « V 2 , F ol (x) « F ,(X) and G tj (k,X) « G\f(k). Then, 
eq. ( |3.47| ) is just the Wigner transform of eq. ( |3.53| ). 
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To conclude this section, let us remark that 5G(k,X) is of the form (compare to 
eq. ([2331) in SQED): 

8G ab (k,X) = p (k)5N ab (k,X) 

= 2n5(k 2 ){6(k )5N ab (k,X) + 9(-k )5N ba (-k,X)}, (3.54) 

where the structure of the second line follows from the first symmetry property ( |3.22|) . 
and the density matrix 8N ab (k, X) satisfies the equation ||] (with t> M = (l,k/fc)): 

dN 

[v ■ D x , 5N(k, x)] = -gv E(x) — , (3.55) 

ak 

which may be seen as the non-Abelian generalization of the Vlasov equation. Note also 
that eqs. ( |3.52|) and ( |3.54j ) hold for both 5G < and 5G > , which are equal in the mean field 
approximation: 

5G < (k,X) « 8G > (k,X) = 5G{k,X) = Po {k)5N{k,X). (3.56) 

This results from the fact that the spectral density p(k, X) = G y (k, X) — G < (k, X) is not 
modified in the present approximation: p(k,X) « po(k) = 2iie(kQ)5(k 2 ). 



3.5 Collision terms in QCD 

As we have seen in the previous section, the colour background field induces a fluc- 
tuation 5Gij ~ (Dx/T)Gq in the Wigner function of the hard transverse gluons. For 
Dx ~ 9 2 T, this fluctuation is of order g 2 G , and the various terms in eq. ( |3.50|) are all of 



order D 2 x Gq ~ g 4 T 2 Go. In this case, the collision terms cannot be neglected and must be 
added in the r.h.s. of eq. ( |3.50| ). 



In order to compute these terms, we consider, as usual, the difference of the self- 
energy terms in the r.h.s. of eqs. ( |3.19| ) and ( |3.2U| ). These involve convolutions of self- 



energies and propagators which yield, after a Wigner transform, 

J d 4 zG(x,z)£(z,y) -^G(k,X)V(k,X) + *-{G, S} pi? + (3.57) 

up to terms involving, at least, two soft derivatives. The Poisson bracket {G, £}p_b is 
defined as in eq. ( 2.44 ). To simplify writing, we have left aside the Minkowski indices; 



these will be added on the final equations. Note also that £ and G are colour matrices, 
so their ordering is important. 



Collecting all the terms without soft derivatives in the r.h.s. of eqs. (|3.19|) and 
7201), we obtain: 



C(k,X) = i[^ R G < -G^a + ^Ga-GrE' 

= -i({G > ,S < }-{S > ,G < })-z[ReE fl ,G < ]+«[ReG' i? ,S < ], (3.58) 
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where the various parantheses stand for colour commutators or anticommutators. In 
writing the second line above, we have also used the relations ( p. 46 ). 



We now proceed to some approximations. Recall first that both the soft gradi- 
ents, and the amplitudes of the background fields and of the fluctuations SG or SH, are 
controlled by powers of g. Writing for instance G > (k,X) = G^ q (k) + SG > (k,X) and 
£<(£;, X) = £<(*) + 8E<(k,X), we have SG> ~ g 2 G>, and similarly <5£< ~ # 2 E< (see 
the next section for the latter estimate). Thus, to leading order in g, we can linearize 
C(k, X), eq. (|3.58|) , with respect to the off-equilibrium fluctuations. Since the equilibrium 



two-point functions are diagonal in colour (e.g., = 5 ab G eq ), the two commutator terms 
in eq. ( |3.58| ) simply vanish, while the anticommutator terms yield: 

C(k,X) ~ - (G>8E< + *G>E<) + (SE>G< + E>6G<). (3.59) 



Each of the terms in the above equation is of order g 2 G eq Yi eq . At this order, the Poisson 
bracket in eq. ( |3.57| ) can be neglected. Indeed: 

{G, s} pB = d k G ■ d x Z - d x G ■ d k Z = d k G ■ d x SZ - d x 5G ■ <9 fe £, (3.60) 

where we have used, e.g., d x T,(k,X) = d x ST,(k, X). With <5£ ~ g 2 ^eq and a similar 
estimate for SG, each of the two terms above is ~ g 2 (d x /T)G eq T leg ~ g 4 G eq T, eg . 

Thus, at the order of interest, the only relevant collision terms are those displayed 



in eq. ( 3.59 ). This corresponds to the quasiparticle approximation introduced in Sec. 2.4. 
Indeed, it can be verified that with the Poisson brackets excluded, the hard gluon spectral 
density p(k, X) satisfies the same equation as in the mean field approximation 0], so that 
p(k,X) ~ Po(k) = 27ie(ko)5(k 2 ). This has the consequence discussed at the end of the 
previous section, namely: 

5G<(k,X) ss 5G > {k,X) = 5G{k,X) = 27r5{k 2 )e{k )5N{k,X), (3.61) 



and the density matrix 5N(k, X) has the structure displayed in eq. ( 3.54|) . Here, however 



5N(k, X) will be shown to satisfy a Boltzmann-like equation, with the collision terms in 
eq. (|3.59| ). In the same approximation, the equilibrium 2-point functions G> and Gf q 
coincide with the free respective functions, as given in eqs. ( |3.15| ) and Q2.14 ). 



We end this section by completing the following two tasks: (i) First, we shall rewrite 
the collision terms ( |3.59| ) in a manifestly gauge-covariant way. (m) Then, we shall specify 
the tensor structure of the collision terms in Minkowski space. 

For point (i), it is enough to replace the non-covariant fluctuations 5H and SG in 
eq. (|3.59|) by the corresponding gauge-covariant expressions <5E and SG (cf. eq. (|3.31|) ): 



SG(k,X) = SG{k,X)-g(A(X)-d k )G eq (k), 

SZ(k,X) = S±(k,X)-g(A(X)-d k )'E eq (k). (3.62) 
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This yields: 



C(k,X) 



G>6E< + SG>E< 



SX>G< 



E> q SG< 



(3.63) 



which turns out to be the same expression as above, eq. ( |3.59| ), except for the replacement 
of ordinary by gauge-covariant Wigner functions: The corrective terms in eq. (|3.62| ) do 
not contribute to C(k, X) since they are proportional to the collision term in equilibrium, 
which is zero: 



g(A(X).d k )(G> q E< - G<E> 



eq eq 



0. 



(3.64) 



And, actually, eq. (|3.63|) is formally the same as in SQED (cf. eq. (|2.57|) ) 



Concerning point (m), recall that the equilibrium 2-point functions Gf ^ u and G^ v 
- which coincide here with the corresponding tree-level functions; cf. eq. ([3.15 ) - 
have only spatial, and transverse, components. These will in turn project the tranverse 
components of the self-energy fluctuations Stlfj and <5S^- in the collision term 



Accordingly, the kinetic equation for 5Gij(k,X), which reads (cf. eq. ( |3.50| )) : 

2[k ■ D x , 5G {j (k,X)] ~ 2(% - k i k j )9k a F af3 (X)d /3 G (k) = C^X), (3.65) 
admits a transverse solution: 

8Gij(k, X) = (5ij - kik^SGik, X) , (3.66) 
as in the mean field approximation. Defining transverse projections in the usual way, e.g., 



5t>(k,X) = (1/2) (fy - k i kj)S'Ef j (k,X), 
we are finally led to the following kinetic equation for SG(k, X) : 



k-D x , 8G(k,X)] - 2gk a F Q ^X)d p G (k) = C(k,X), 



with the collision term: 



(3.67) 



(3.6* 



C(k,X) = -T T (k)5G(k,X) + (5S>G< - <5£<G>j. (3.69) 
We have recognized here the equilibrium damping rate for the transverse gluons (cf. 



eq. (|239l) ) 



IY(fc) = (£< (*) - E> (k)) T . 



(3.70) 



(Note that, in what follows, the subscript T on transverse quantities will be often omit- 
ted.) Eqs. ( |3.68| )-( |3.69| ) are manifestly covariant under the gauge transformations of the 
background field. 
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3.6 The hard gluon self-energy out of equilibrium 



In this subsection, we use perturbation theory to compute the transverse gluon self-energy 
5Ti(k, X) to the order of interest. Specifically, we shall find that ~ g 2 T, eq ~ g A T 2 , so 
that the collision terms in the r.h.s. of eq. (|3.68 ) are of the same order as the drift and 
mean field terms in the l.h.s. 

We start with the ordinary (i.e., non-gauge-covariant) self-energy H(x,y) out of 
equilibrium. As in SQED, the leading-order collision term corresponds to scattering via 
one gluon exchange, as illustrated in Fig. |I[ However, as stated in the Introduction, 
we are mostly interested here in colour relaxation, for which the relevant collisions are 
dominated by soft momentum transfers, g 2 T ^ q ^ gT (cf. Sec. 3.8 below). Accordingly, 
the virtual gluon in Fig. |1] is always soft and the self-energy which describes this collision 
is the effective one-loop diagram depicted in Fig. f|. This is formally the same diagram 
as in SQED, except that, now, the continuous line in Fig. ^ refers to a hard transverse 
gluon and the wavy line to a soft virtual one (which can be longitudinal or transverse). 
The bubble on the wavy line denotes, as usual, the resummation of one-loop polarization 
tensor in the propagator of the soft gluon. 

In thermal equilibrium, the hard line in Fig. |] is a free propagator, while the soft one 
is the HTL-resummed propagator, as introduced in Sec. 2.6. (In the HTL approximation, 
the gluon and photon propagators are formally the same up to the replacement of the 
Abelian Debye mass m 2 D = e 2 T 2 /3 by the non-Abelian one, m 2 D = g 2 N c T 2 /3 0,0.) Then, 
the self-energy in Fig. (| yields a contribution of 0(g 2 T 2 ) to the thermal interaction rate 
T(fc), eq. ( 3.70|) . Thus, the first collision term in eq. (|3.69|) can be estimated as: 



T(k)5G(k,X) ~ g 2 T 2 5G ~ g 4 T 2 G . (3.71) 



Consider now the other terms in eq. ( [3.69D , which involve the off-equilibrium self-energy 



ST,. A typical term contributing to SH is obtained by replacing the equilibrium propagator 
Go m the hard line in Fig. [| by the respective off-equilibrium fluctuation SG. Thus, 

5£ ~ (5G/G )Z eq ~ g 2 T, eq ~ g 4 T 2 , (3.72) 



which contributes to the collision terms at the same order as the damping rate in eq. (|3.71|) . 



Moreover, the off-equilibrium effects enter also the soft gluon line in Fig. |], via the po- 
larization tensor. This will be computed in the next section, where we shall see that the 
net effect is also of order g 4 , as in eq. (|3.72j ). 



Let us turn now to the explicit evaluation of E. To this aim, we need the three-gluon 
vertex between two hard gluons and a soft one in the presence of the background field. This 
can be read on the Yang- Mills action, in the following way: Split the fluctuating gluon field 
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(that we have originally denoted as a%(x); cf. eq. ( |3.5| )) into soft and hard components, 
and use different notations for the two. That is, replace^ a%(x) — > a%(x) +A%(x), where 
the new field a p is hard (it carries momenta k ~ T), while the field A p is soft (with typical 
momenta g 2 T C g < T). There is no difficulty with this separation (e.g., no problems 
with gauge symmetry) since the collision terms will be saturated by soft momenta q ^ gT. 
The Yang-Mills piece of the action reads then: 

Sym = -J c d*x^(F; u [A + a + A]) 2 

= ^J c d 4 xaZ{g IMU D 2 -D^D u + 2igF IMU } ab a» + ... (3.73) 

where = D^[A + A] = -D M [v4] + igA^, and the dots stands for terms of cubic or 
quartic order in a M , which are unimportant here. We still have to isolate the trilinear 
couplings a^a u A p from the equation above. After some algebra and integration by parts, 
the relevant interaction piece of the action is obtained as: 

Sj = ig [ d A xa»Jg^(A-D) - \{A^D V + A V D^) + 2[D„, aA <, (3.74) 

JC I I ) ab 

where now = -D M L4] = d^ + igA^ is the covariant derivative defined by the background 
field alone, and all the fluctuating fields are explicit. 



Only the first two terms in eq. (|3.74 ) will be important. In these terms, the covariant 



derivatives act on the hard fields a p and give rise to vertices with hard momenta. The 
third term, on the other hand, involves a covariant derivative acting on the soft field A 1 , 
and is sub leading (by a factor of q/k ^ g). In what follows, we shall ignore this term, 
and focus on the self-energy built out of the first two terms in Sj. We write: 

Si = |r^ pA J c d 4 xa»A p D x [Ay, (3.75) 

with an implicit trace over the colour indices (the symbol T^ up x has been defined in 
eq. ( 3.18| )). With the three-particle vertex above, it is a straightforward exercise to con- 



struct the self-energy displayed in Fig. ffl. This reads: 

E f v (x, y) = - 4 x t. rV^r^ (T a ) cd (T b )- C(I V^(x, y) (d p G^(x, y)D}%, (3.76) 
where G a/3 (x,y) is the hard gluon propagator and V lS (x,y) the soft gluon propagator: 

V$(x,y) = (T c Al(x)At(y)). (3.77) 

Finally, the factor 4 takes into account the fact that, strictly speaking, there are three 
other terms similar to the one above, which yield the same contribution to the order of 
interest (see below). 

■^Note that we preserve the notation A£(x) for the colour background field. 
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By choosing the time variables in eq. ( |3.76| ) on opposite sides of the contour, we 



deduce expressions for the Wigner functions and E < . The upperscripts > and < will 
be often omitted, to simplify writing. 

Consider first the equilibrium limit of eq. ( |3.76|) , where = and the internal 
propagators are unit matrices in colour: By using {T a ) c d (T b ) cc i = —N c 5 ab and going to 
momentum space, we obtain (E eg )^(A;) = 5 a bT,^(k) with : 

= g 2 N c T^ p T uP&x J-0^( k ~ ^ ~ V) XV 2 S M G ^( k ~ ?)■ ( 3 - 78 ) 

In the construction of the collision terms, we shall need only the difference £ < — S > (cf. 
eqs. ( |3.69| ) and ( |3.70| )), and the resulting integral will be dominated by soft momenta. 



With this in mind, we shall neglect q next to k ~ T in the vertices in eq. ( |3.78| ). Further- 
more, to the same order, G eq ~ Go (which has only spatial and transverse components; cf. 
eq. ( p,15| )) and V eq « *T> (which is the HTL-resummed gluon propagator; cf. Sec. 2.6). 
We finally get the following estimate for the transverse gluon self-energy in equilibrium: 

E<(*) = (l/2)(5 ij -k i k j )^ < = Ag 2 N c k"k x J *V< x (q)G<(k-q), (3.79) 

together with a similar expression for S^(fc). The following identity has been useful in 
performing the Minkowski algebra (with Vijijt) = 5ij — kihj): 

(l/2)7> -(k) T lhp kf>V lm (k) Tj m $\k x = Ak y k s . (3.80) 



We now return to the general expression in eq. Q3.76P and evaluate the off-equilibrium 



fluctuation 5S. Since we consider only small deviations away from equilibrium, we can 
linearize this expression, as we have already done for the collision terms in Sec. 3.5. We 
thus get, keeping explicit only the colour indices: 

5E ab (x,y) = (T a ) cd (T b ) 5j {5V dd (x,y)(d^G (x-y))6 cS 

+ 6 d /V(x - y) 5(DP x G(x,y)D^)J. (3.81) 

From this, we shall construct the gauge-covariant self-energy <5£ a6 (s, X), as explained in 
Sec. 3.3. First, we replace the non-covariant fluctuations in the internal propagators 
8T> and SG in terms of the corresponding gauge-covariant fluctuations 5T> and SG (cf. 
eq. (|3T3q ) and 



6V{x,y) = 6V(s,X)-ig(s.A(Xj)*V(s), 
5{D^y)D X y ] ) = -d&5G(s,X) + ig(s-A(X))d&G . (3.82) 

Then, we define the covariant self-energy as in eq. ( |3.62|) : 



S±(s,X) = 5H{x,y) + ig(s-A(X))X e Js), (3.83) 
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with S e g of eq. ( |3.78| ). The result of these operations has the rather simple form: 

5± a \s,X) = (T a T b ) cd {6V cd (s,X)d&G (s) + *V(s)d&6G cd (s, X)} (3.84) 
so that, after a Wigner transform: 



S&*(k,X) = -(T a T b ) cd J -0^(k - q)»(k - q) : 



x{6V cd (q,X)G (k -q) + *V(q)5G cd (k - q,X)}. (3.85) 
By putting back the factors of g 2 and the Minkowski indices, we finally obtain the following 



expression for the gauge-invariant self-energy fluctuation 5t^ v : 



[IV 



5±; b u (k, X) = g 2 r w r^ A (T a T b ) cd jS^{k- q) p (k - q) X 

x {SVS(q, X)Gf(k -q) + *V^(q)5G a J(k - q, X)}. 

(3.86) 

This expression is very close to the corresponding expression in equilibrium: eqs. ( |3.78j ) 
and (|3.86|) involve the same momentum-dependent vertices, and the equilibrium propa- 



gators of eq. ( [3.781 ) have been simply replaced in eq Q3.86|) by the (linearized version) of 



the respective off-equilibrium propagators. The only significant difference is the colour 
structure, which is trivial in equilibrium. 

As before, we can neglect the soft momentum q in the vertices of eq. (|3.86| ), and 
take the transverse projection of this expression, to obtain: 

8±$(k,X) = 4g%k x (T a T b ) cd J ^{*V^q)5GUk-q,X) + 5V^(q,X)G (k-q)}, 

(3.87) 



which is the non-equilibrium generalization of eq. ( [3.79 ). 



To conclude this section, let us return to a previous remark according to which the 
complete self-energy S should involve three other terms in addition to the one displayed 
in eq. ( |3.76| ). In these terms, the covariant derivatives act differently on the two internal 
propagators, which then results in modification of the momentum-dependent vertices. 
For instance, we meet terms like D%Dy [G(x, y)D(x, yu which, after covariantization and 
Wigner transform, yield the same result as in eq. (|3.86|) , except for the replacement 
(k — q) p (k — q) x — > k p k x in the vertices. However, this difference is not important here 
since we neglect q next to k in the vertices. The same holds for the other two terms, so 
that the total contribution is, indeed, four times the contribution of the term displayed 
in eq. ( |3.76|) . Hence the factor 4 in eq. Q3.76|) . 
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3.7 The off-equilibrium propagator of the soft gluon 



The above expression for 5T,(k,X) involves the off-equilibrium propagator of the soft 
gluon 6V pX (q,X), to which we turn now. The relevant off-equilibrium effects are encoded 
in the soft gluon polarization tensor, which we denote by Tl a £ v (x,y). In equilibrium, this 
reduces to the corresponding hard thermal loop 8 Tl e ^{x — y). Thus, the calculation 
below provides a generalization of the HTL polarization tensor out of equilibrium. 

The Kadanoff-Baym equations for the soft gluon propagator T> lJih ,{x,y) are formally 
identical to those for G^(x,y), i.e., eqs. ( |3.19|) and ( |3.20| ). For instance: 

(g p D 2 - D /Ji D p + 2igF f f) ^ V< - J d 4 z g, x II jf V< = J d 4 z g Xu H< p V A \ (3.88) 

We shall also need the retarded propagator V R u (x,y), which obeys (cf. eq.( |2.24j )): 

(g p p D 2 - D P D P + 2igF§V R v {x,y) - J d 4 zg Xp Il R x V% = g^5^(x - y) . (3.89) 

A priori, the self-energy U pu in these equations involves both interactions with the hard 
fields a p and self-interactions of the fields A p . However IL^ will be dominated by the 
one-loop diagrams depicted in Fig. ^ where the internal lines are hard. 

As in the Abelian case, the two equations above imply a relation between T> < and 
n< (cf. eq. ( P^D ) : 

^paO^ V) = - J d A z 1 d 4 z 2 g pa V% l (x, z x ) U<,(zi, z 2 ) V v /(z 2 , y)g px . (3.90) 

A similar relation holds in between D > and Il > . In particular, in thermal equilibrium, 

*V<M = - (*V R (q) H< (q) *V A (q)) pX , (3.91) 

where *T>(q) is the HTL-resummed propagator. 

Out of equlibrium, we can compute the Wigner transform of eq. (|3.90| ) to get: 

V< x (q,X) = - (V R (q,X)U<(q,X)V A (q,X)) x . (3.92) 



This equation holds up to corrections of 0(dx/q)- This is an important limitation here 
since dx ~ g 2 T, while g 2 T ^ q ^ gT. However, since it will turn out that the colour 
relaxation rate is only logarithmically sensitive to the ultrasoft momenta q ~ g 2 T, it can 
be argued that the terms which have been neglected in the above gradient expansion are 
suppressed by a factor of 1/ ln(l/g). In this sense, eq. ( |3.92|) is still correct to logarithmic 
accuracy. 

The difficulty we are facing here comes from the necessity to perform a gradient 
expansion in the presence of long range interactions, and has already been alluded to 
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in the case of SQED in Sec. 2.5 (see the discussion after eq. ( j2.7C| )). In principle, one 



could develop a more accurate approximation scheme by allowing for a collision term 
non-local in X. Specifically, in the exact equation ( |3.90| ), we can safely treat x and y as 
neighbouring points, with \x— y\ ~ since the propagator T> < (x, y) is to enter the hard 
gluon self-energy y) in eq. (|3.76|) ; we can proceed similarly with the points Z\ and z 2 : 
\z\ — z 2 \ ~ 1/T, since the polarization tensor H<(zi,Z2) is dominated by the hard loop 
(e.g., p is a hard momentum in eq. ( |3.97|) below). But the points x and z\ (or y and z%) are 
relatively distant one from the other, since they are related by the long-range propagator 
T>r(x, zi) (respectively, by T>a{z 2) y))- Thus, a more accurate gradient expansion should 
treat X = (x+y)/2 and X' = (z 1 + Z2)/2 as distinct points (the "end-points" of the virtual 
gluon line in Fig. |l|), which would then lead to a collision term which is non-local in X'. 
However, the construction of such a non local collision term goes beyond our present goal, 
and we shall stick to the local expression in eq. ( |3.92|) . 

The gauge-covariant fluctuation 5T> < (q, X) is obtained from eq. (|3.92|) by first 
linearizing with respect to the off-equilibrium fluctuations (e.g., Il < (g, X) = Uf q (q) + 
8H < (q, X), with 5U < ~ g 2 Uf q ), and then replacing the ordinary 2-point functions with 
the corresponding covariant ones. This gives: 

6V<(q,X) = SV<(q,X) + g(A(X) ■ d q )*V<(q) = -*V R (q) 5U<(q, X)*V A (q) ...,(3.93) 

where the dots stand for terms involving the off-equilibrium deviations of the retarded, 
or advanced, functions, but the equilibrium self-energy Hf q (q). It is easy to verify that 
such terms will eventually cancel in the collision terms (as they do in the Abelian case: 
cf. the remark after eq. ( |2.80| )), so we shall ignore them in what follows. 

We thus need the off-equilibrium polarization tensor 511 < (q, X) for soft q. This is 
determined by one-loop diagrams which look formally as in SQED (cf. Fig. |6|), except 
that, in QCD, the internal lines denote hard transverse gluons. The tadpole diagram in 
Fig. |[b does not contribute to the collisional self-energies Il > and Il < , but only to the 
retarded self-energy Il#, which we know already to be the HTL (recall that we only need 
this in equilibrium; cf. eq. ( |3.93| )). Therefore, in what follows we shall focus on the non- 
local^ self-energy in Fig. |].a, which we evaluate by using the three-particle vertex a^a u A p 
in eq. ( |3.75|) . This yields: 



V) = \ (T a ) cd (T%-{ (D^(x, y))JcP{x, y)D}1) a 

xG%(x,y)(D^ 5 (x,y)D^)J. (3.94) 



g ln Coulomb's gauge, there is another tadpole coming from the diagram in Fig. |^.a where one of the 
internal lines is transverse, and the other one is longitudinal (and static) p8[ . This too contributes to 
the retarded/advanced propagators, but not to the collisional self-energies. 
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Starting with this expression, the linearization and the covariantization proceed along the 
same lines as for the hard gluon self-energy in eq. (|3.76|) . (This is legitimate since the loop 
integral is dominated by hard momenta, so that U^x, y) is localized at \x— y\ ^ In 
this process, we use identities like eq. ( |3.30|) , ( |3.32| ) and (|3.82| ) to replace ordinary Wigner 
functions by gauge-covariant ones, and define the covariant polarization tensor as usual: 

5tl(s,X) = 5U(x,y) + ig(s-A(xj)ll eq (s), (3.95) 

with 

n» = ^ T, paP T vXlS J j0p P P P X G^(p + q) Cf (p). (3.96) 

Here again, we have neglected q ^ gT next to p ~ T in the vertices. 

The final result is quite predictible: The covariantized fluctuation Stl(q, X) is for- 
mally similar to the equilibrium self-energy ( |3.96j ), except for the replacement of the 
equilibrium internal lines by (linearized) off-equilibrium gauge-covariant Wigner func- 
tions. Once again, this simple result holds only after covariantization, and implies that 
the internal momenta in eq. ( |3.97|) below should be interpreted as kinetic momenta (recall 
the discussion after eq. (|3.86|) ). Specifically: 

SfftfaX) = ^T, pa ,T^ /5 (T a T b ) cd J '-^P P P X 

x{G^(p + q)6G s f c (p,X) +5G a c ](p + q,X)Gf(p)}. (3.97) 

At this point, we remember that the hard gluon transverse functions (in or out of equi- 
librium) are purely spatial and transverse, so that the above equation can be further 
simplified to: 

(StintMX) = ^9 2 (T a T b ) cd J-0^p^ 

x{G>(p + q)5G< c (p,X) + 5G> d (p + q,X)G<(p)}, (3.98) 
where we have reestablished the upper scripts > and <. 



3.8 The Boltzmann equation for colour 



We are now in position to explicitly compute the collision terms in eq. ( |3.69|) . We write 
C = Ci + C 2 , with 

Ci = —F(k)8G(k,X), C 2 = 5± > {k,X)G<{k) - 5± < {k,X)G>{k). (3.99) 
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The first piece C\ involves the (equilibrium) damping rate for hard transverse gluons T(k) 
which can be computed from eqs. Q3.79Q and (|3.91|) above: 

k?k\*V>M G>(k -q)- *D< (?) G<(k - q) 

d 4 q 



r(k) 



4N c g' 



(2tt)< 



k p *V%(q) *Vf(q)k x \ (II> (g)G£(* - q) - U< u (q)G^(k - q)). 

(3.100) 



Here, II 



11^ is the polarization tensor in equilibrium, as given in eq. ( |3.96| ). By 



inserting it in eq. ( |3.100|) , one finds (with p' = p + q and k' = k — q): 

T(k) = Nl [ \M\ 2 {G<(p)G>(p')G>(k') - G>(p)G<(p')G<(k>)}, (3.101) 



where 



d 4 q f d 4 p 



pg J (27T) 4 J {2TTP 

and we have recognized the matrix element squared for the collision depicted in Fig. |I] : 



\M\ 2 = 16/ k;Vf(q)p, p v *Vf{q)k 



(3.103) 



The second piece C 2 in eq. (|3.99| ) can be similarly computed by using eqs. (|3.87|) , (|3.93| ) 
and fl3.98j ). One gets: 

Ct = -N c (T a T b ) cd [ ^^^(^^(pOg^w-^c^g^^g^^)}^^^, 



Cab 
99 = 



PI 



G 



(T a T b ) c - c (T c T 5 U f \M\ 2 {G<(k)G>(k')G>(p') - G>{k)G<{k')G<{p')}5G dd {p,X), 

Jpq 

\M\ 2 {G<(k)G>(k')G<(p) - G>(k)G<(k')G>(p)}5G^(p',X). 

(3.104) 



23 = -(T a T b ) c5 (T c T c 



In writing these equations, we have used the fact that £>G < ~ 5G > = 5G in the present 
approximation (cf. eq. ( |3.61|) ). The piece C21 comes from the first term in the r.h.s. 
of eq. ( |3.87| ), which describes fluctuations in the hard propagator inside H(k,X) (the 
lower line in Fig. ^). The other two pieces, C22 and C23, come from the second term in 
eq. (|3.871) an d describe fluctuations in the soft (upper) line in Fig. |]. Clearly, these three 
terms C21 C22 and C23 are associated with fluctuations along the external lines "to be 
summed over" in Fig. |] — namely, the lines with momenta p, p' and k' — , as opposed to 
C\ which describes fluctuations along the incoming line with momentum k. 

We shall verify in a moment that the phase-space integrals in eqs. ( |3.101p and 
( |3.1041 ) are indeed dominated by soft exchanged momenta q, which justifies our previous 
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approximations. This allows us to make some further simplifications, as follows: Recall 
first that all the Wigner functions in these equations are distributions with support on 
the tree-level mass-shell (cf. eq. ( ft.61| )). E.g., 

G<(k')=p (k-q)N(k -q ), 5G(p',X) = p (p + q)5N(p + q,X), (3.105) 

with po(k) = 27re(/co)<5(fc 2 ). In such expressions, we cannot neglect q next to k (or p) 
within the mass-shell 5-functions, but we can still do that in the occupation numbers: 

G<{k')~Po(k-q)N(k ), 5G(p',X)^p (p + q)6N(p,X). (3.106) 

This yields, e.g., (compare to the Abelian expression in eq. ( 2.82| )) : 



T{k) ~ N 2 j dT\M\ 2 {N(p )[l + N(p )][l + N(k )}~[l + N(p )}N(p )N(k )} 

= N 2 C JdT\M\ 2 N( Po )[l + N(p )], (3.107) 

and the phase-space measure / dT has been defined in eq. ( |2.80[ ) . Similarly, 

C£ ~ p Q (k)N c (T a T b ) cd 6N cd (k,X)jdT\M\ 2 N( Po )[l + N(p )}, (3.108) 

which is of the same form as the damping rate contribution C% = — T{k)8G{k, X) (cf. 
eq. ( 3.1071) ), and can be combined with the latter to yield (in matrix notations): 



C l + C 2l = - Po (k)^[T a ,[T a ,5N(k,X)]] j dT\M\ 2 N(p )[l + N( Po )]. (3.109) 
The remaining two terms C22 and C23 can be similarly simplified to yield: 
C? 2 + C£ ~ p (k)(T a T b ) c - c (T c T- c ) dd :J dT \M\ 2 N(k ) [1 + N(k )] (6N dd - 6N dd ) (p, X) 

AT „ 

= Po (k)^(T% b j dT\M\ 2 N(k )[l + N(k )]Tr{T e 5N(p,X)), (3.110) 

where the second line follows after some elementary colour algebra. 

Thus, unlike the Abelian case in Sec. 2.6, where the IR contributions to the four 
terms in eq. ( |2.79| ) cancel each other, here we do not have a complete cancellation because 
of the non-trivial colour structure. By putting together eqs. ( |3.109| ) and (|3.110|) , we finally 
derive the following expression for the non- Abelian collision term: 

C[5N] = -p (k)^JdT\M\ 2 {N( Po )[l + N( Po )}[T a ,[T a ,5N(k,X)]] - 

-N(k )[l + N(k )]T a Tr(T a 5N(p,X))}. (3.111) 



Together, eqs. ( |3.68| ) and ( |3.111| ) determine a Boltzmann equation which describes colour 



relaxation in hot QCD. By also using eq. Q3.54Q , this can be rewritten as an equation 
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for the density matrix 8N a b(h, X). To this aim, we need the positive-energy projection 
(ko = Ek = |k|) of the collision term. For soft q and k = Sf, ~ T, k' is positive as well: 



Po(k')\ 



2tt 



fc =£fc 



2s 



k — q 



8{s k -q0~ £k~q) 



2tt 
2s~ k 



5{q - q- v) 



(3.112) 



where v = k is the velocity of the incoming particle with momentum k. Concerning po, 
this can be either positive, po = e p , or negative, po = —s p , and the two situations yield 
identical contributions. We thus replace (with v' = p): 



Po(p)po(p + q) 



( 2n \ 2 

2 I 2^-1 8(p -£ p )8{q -q-v'). (3.113) 



Also, for on-shell momenta ko = £k and po = £ p , the matrix element (|3.103|) becomes: 

|A^| 2 = Wg'elel + (q x v) • (q x v') *V t (q)\\ (3.114) 

with the same notations as in the Abelian case (cf. eq. ( |3.114| )). 

By also using the identity N(p)[l + N(p)] = —T(dN/dp), we are finally led to the 
following Boltzmann equation, which is the main result in this paper: 



[vD x , 5N(k,X)} + 9 vE(I)^ = C[6N], 



(3.115) 



with the collision termQ: 

d 3 p 



C[6N}=g*N c T J -0^ $(v ■ v'){^ [T<\ [T a , 5N(k, X)}} - ^rTr(HiV(p,I)) 



(3.116) 



and the collision integral (v t and v£ are the velocity projections transverse to q): 

$(v • v') = (2tt) 2 J ^0 5(q - q ■ v)<% - Q ■ v')|*^(?) + (v* ■ v' t ) *V t (q)f (3-117) 

The collision terms above are identical to those written down in Ref. on an heuristic 
basis (cf. eq. (3.26) in Ref. 0). 

The above equations can be further simplified by noticing that the corresponding 
solution 8N(\a, X) can be written in the form: 



8N(k,X) = -gW(X,v) (diV/dife), 



(3.118) 



where the new function W(X, v) (a colour matrix) depends upon the velocity v (a unit 
vector), but not on the magnitude k = |k| of the momentum. This function satisfies a 



h Note that the overall normalizations in eqs. (3.111) and (3.116) arc different. 
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simpler equation (we change X — > x from now on since this is the only space-time variable 
left in all the equations to come): 

a 2 T r d.0' 

[v-D x ,W(x,v)) = v-E(x) - m^^-y — $(v- v') 

x{[T a , [T a ,W(x,v)}} - T a Tr (T a W(x,V)}, (3.119) 

where the angular integral runs over all the directions of the unit vector v', and the Debye 
mass m 2 D comes out after performing the radial integral over p = |p| : 

2 9 2 N C roo 2 dN g 2 N c T 2 
m D = — / dpp z — = . (3.120) 



7r 2 jo dp 

The factorized structure of the colour density matrix in eq. ( 3.118| ) has been first recog- 
nized at the level of the mean field (or Vlasov) approximation (cf. eq. ( |3.55|) ), where it 
has been shown to have a simple physical interpretation ||. It is remarkable that such a 
structure persists after the inclusion of the collision terms. 

Consider furthermore the colour structure of eq. ( |3.119| ). In constructing the in- 
duced current ( [3.1|) , we only need the first colour moment of the density matrix W(x, v), 
namely W a (x, v) = (l/iV c )Tr (T a W(x, v)). The equation satisfied by W a (x, v) follows 
from eq. fl3.119| ) by taking the appropriate colour trace, and reads: 

(v ■ D x ) ab W b (x,v) = vE«W-m^ J ^®(v -V){w a (x,v) -W a (xy)}. 

(3.121) 

This is the equation which has been announced in the Introduction (cf. eq. ( |1.2|) ). The 
collision term in this equation has a simple physical interpretation. We may indeed 
identify 



m 2 D 



g 2 N c T fdn'. 9 2 N c Ts m D 



*(v.v<)= 7 =<4^(ln^ + 0(l)), (3.122) 

47T AlX V /i ' 

which is the damping rate 7 = T(k = k)/4k for a hard transverse gluon with velocity 



v 



3lj |32[| . In the present approximation, the gluon damping rate suffers from the same 
logarithmic IR divergence as the damping rate for a charged particle in the Abelian case 
(cf. Sec. 2.6). In the equation above, this has been cutoff by hand, by introducing an IR 
cutoff fi (see the r.h.s. of eq. fl3.122j )). 

In previous studies of the colour conductivity or the damping rates, it has been 
generally assumed that the IR cutoff is provided by non-perturbative magnetic screening 
at the scale g 2 T 0, ||. This yields a damping rate 7 ~ aN c T ln(l/g) (with a = g 2 /47i), 
but the constant term under the logarithm cannot be determined: indeed, this is sensi- 
tive to the details of the magnetic screening, which remains poorly understood. In the 
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framework of the effective theory recently proposed by Bodeker ||, \x should be rather 
understood as an intermediate scale /i ~ g 2 Tln(l/g) separating the perturbative physics 
at "hard" (k ~ T) and "semi-hard" (/i ^ q ^ gT) momenta from the non-perturbative 
physics at "ultrasoft" (q ^ g 2 T < fi) momenta. Then, [i will hopefully cancel in any 
complete calculation, via a matching between the perturbative and the non-perturbative 
(e.g., lattice) calculations. However, our above derivation of the collision term shows that, 
within the present formalism, it is not possible to go beyond leading-log accuracy (which 
makes the issue of matching superfluous): this reflects the present limits of the gradient 
expansion in the presence of long range interactions (cf. the discussion after eq. ( |3.92|) ). 



Within this logarithmic accuracy, the integral in eq. ( |3.117| ) can be performed by 
replacing the matrix element for transverse scattering by its infrared limit as go ^ <? ~~ *• 0, 
namely (see eq. ( |2.84j )) 



rV t (q ,q)\ 2 ~ \ — 4r^M. (3 .i 2 3) 

One thus obtains: 

2 (v ■ v') 2 

$(v • v') ~ 1 ] ln(l/</) , (3.124) 

TT Ttljj . I _ V V' 2 



where the ln(l/g) comes from (compare to eq. ( p.85|) ) : 



mD da aT 1 

^ ~ In y — ~ In -, (3.125) 

Q V 9 



with n ~ g 2 T\n(l/g). Then, eq. (|3.121|) reduces to (with 7 ~ aN c T\ia(l/ g)) : 



[ 4 r dQ' (v ■ v') 2 I 

[v ■ D x )W(x, v) = v E(x) - 7 <^ W(x, v) - - / . v ; W(x, v') } .(3.126) 

71 J An . 1 _( v . v >)2 



This is precisely the kinetic equation which generates Bodeker's effective theory || 0. 

The previous equations have been obtained by working in the Coulomb gauge for 
the hard field fluctuations (cf. eq. (|3.12|) ), but we expect the final results — namely, 
the Boltzmann equation in eqs. ( |3.115| )-( |3.116| ), or (|3.121|) — to be actually gauge-fixing 



independent. Except for the collision term, this has been explicitly proven in 0] (see also 
Refs. f28| , |29|j). Moreover, on physical grounds, the collision term should be gauge-fixing 
independent as well, since it involves only the off-equilibrium fluctuations of the (hard) 
transverse gluons, together with the (gauge- independent) matrix element squared ( |3.114|) . 



(See also the discussion after eq. |3.11| .) 



Furthermore, the previous equations are manifestly covariant under the gauge trans- 
formations of the background fields. This ensures that the induced current (|3.1|), which 
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can be expressed in terms of W a (x, v) as follows: 

= m 2 D J^v»W a (x,v), (3.127) 

transforms indeed as a colour vector in the adjoint representation. Moreover, this current 
is covariantly conserved, 

D,r = 0, (3.128) 



as necessary for the consistency of the mean field equations of motion Q3.16 ). To verify 



eq. ( |3.128j ), use eq. ( |3.121| ) together with the fact that the collision terms vanishes after 
angular averaging: 

r HO r HO' 

J — J —<S>(w'){w a (x,v)-W a ( X y)} = 0. (3.129) 

The quantities W a (x,v) may be regarded as functionals of the mean fields A%(x), 
as given implicitly by the Boltzmann equation (3.121). The current (3.127) itself acts as 



a generating functional for the one-particle-irreducible amplitudes of the ultrasoft colour 
fields. We can formally write (see, e.g., Refs. |], §, §]) : 



j; = U, U A: + \ T%AIA? C + ... (3.130) 

where n M „(P) is the polarization tensor for the ultrasoft [P ~ g 2 T) gluons, is a 
correction to the 3-gluon vertex, etc. For external momenta of order g 2 T or less, the 
amplitudes in eq. ( |3.130| ) are of the same order of magnitude as the "hard thermal loops" 
f28| , p9| , ||, which they generalize by including the effects of the collisions among the 
hard particles. Like the HTL's, the above amplitudes are gauge-fixing independent, and 
satisfy simple Ward identities which follow from the conservation law ( |3.128| ) by successive 
differentiations with respect to the fields A%. For instance: 

p"rv(p) = o, 
PiT^ P (pi,P2,P3) = n„ p (p 3 ) - n W p(p 2 ) . (3.i3i) 

Such identities express the fact that the effective theory at the scale g 2 T, as obtained 
from the Boltzmann equation |J, is gauge invariant. However, unlike the HTL's — which 



in terms of Feynman graphs correspond to one-loop diagrams |28|, |2£| — , each of the 
ultrasoft amplitudes in eq. ( |3.130|) receives contributions from an infinite series of multi- 
loop Feynman graphs, which all contribute at the same order in g. This will be explained 



in the next section and, in more detail, in a further publication |47 . 
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4 Diagrammatic interpretation of the collision terms 



In this section we provide a diagrammatic interpretation of the collision terms. For the 
case of a scalar field theory, such an interpretation has been worked out in much detail 
in Refs. |33|, [34|], where the Boltzmann equation has been actually derived by resumming 
appropriate classes of Feynman graphs. Here, where the Boltzmann equation has been 
constructed by using different technics, it is still interesting to understand what are the 
diagrams which have been effectively resummed in this construction. To this aim, it 
is convenient to go back to eqs. (|3.99p , ( |3.101|) and (|3.104|) where the off-equilibrium 
fluctuations (like, e.g., 5G(k,X)) are unambiguously associated with each of the colliding 
fields in Fig. |TJ 

Consider first the piece C\ = —T(k)5G(k, X) which involves the interaction rate 
r(/c) for the incoming particle with momentum k (cf. eq. ( |3.101| )). This can be moved 
into the l.h.s. of the kinetic equation and combined with the drift term k ■ Dx to yield: 

(2(k-D x ) + T(k))5G(k,x) = 2g(k ■ F ■ d k )G (k) + C 21 + C 22 + C 23 , (4.1) 

with C 2i as given in eqs. fl3.104| ). After the inclusion of T, the drift operator describes 
also the decay of the incoming particle. For instance, the corresponding Green's function 
satisfies (with ko = k and 7 = T/4k; cf. eq. ( |3.122| )) : 

— i(v ■ D x + 2 7 ) A(x, y- v) = 6®(x - y), (4.2) 

with the retarded solution (below, t = x° — y°): 

A R (x,y;v) = iO{t) 6® (x - y - vt) e" 2 ^ U(x, y), (4.3) 

which shows exponential attenuation in time with a rate equal to 27 (compare to eq. ( p. 60] )). 

The inclusion of the interaction rate T(k) in the drift operator amounts to an approx- 
imate self-energy resummation in the hard gluon propagator. At the level of the original 
Kadanoff-Baym equations ( |3.19| ) and ( |3.20|) , this amounts to moving all the terms in- 
volving G K (x,y) (including the terms TirG < and T,aG k ) into the left hand sides of these 
equations. What we would like to argue now is that the other collision terms in the r.h.s. 
of eq. (|4.1|) can be similarly recognized as vertex corrections. 



For more clarity, consider the response of the plasma to an arbitrarily weak electric 
mean field E . Then, non-linear effects are negligible and the relevant response function 
is the polarization tensor: j£(P) — (P) A^(P) , where P ~ g 2 T is the (ultrasoft) 
momentum carried by the background field. In the mean field approximation, and to 
linear order in the electric field, the induced fluctuation reads (cf. eq. (|3.55|)): 



mktP) = ZHI^m*", (4.4) 
v ' v ■ P + it] dk 
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A(P) 




(a) (b) 



Figure 8: Diagrammatic representation of eq. ( |4.5| ): (a) the first-order fluctuation 
SN^(k, P), as induced by the mean field A^(P); (b) the corresponding contribution 
to the polarization tensor rP"(P). The continuous line with a blob denotes the hard 
gluon propagator after the inclusion of the damping rate. 

where the small imaginary part it] with 77 — > + stands for retarded boundary conditions. 

After adding the collision terms, we end up with an integro-differential equation for 
5N (see, e.g., eqs. ( |3.121| ) or ( }4.1| )), whose resolution is non-trivial already in the weak 
field limit. However, in order to see what are the relevant graphs, it is sufficient to consider 
the perturbative solution obtained by iterations. First, it is straightforward to resum the 
equilibrium damping rate T (i.e., the collision term C\): according to eqs. (^4.1|) and (\i.2\) , 
this amounts to including the damping rate in eq. ( [4.4|) . This provides us with the Oth 
order iteration for 8N: 

w <°> (t ,P ) = -«;- E w ^. ,4.5) 

This can be given the diagrammatic representation in Fig. ||a where the blob on the 
continuous line denotes the resummation of the damping rate in the hard gluon propaga- 
tor, and the wavy line with a bubble attached to it represents a mean field insertion. In 
this approximation, the polarization tensor II AtI/ (P) is given be the one-loop diagram in 
Fig. |.b. 

With 5N^ as above, we can compute the first iteration for the collision terms 
C*2i, i = 1, 2, 3. We shall not write down the corresponding formulae, but simply look for 
their interpretation in terms of Feynman graphs. The term C21 in eq. ( p.l04j ) is associated 
with a fluctuation SN(k',X) in the hard, lower, propagator in Fig. |j; its first iteration is 
obtained by using the approximation (|4.5|) for this fluctuation, and has the diagrammatic 
representation in Fig. |9|.a. To the same order, C22 is represented in Fig. ||.b, while C23 
has a similar representation. With these approximations for the collision terms, one can 
compute the first iteration 5N^(k,P), as well as its contribution to U^ U (P): the latter 



is shown in Figs. ff(J.a and pllb. Clearly, Fig. 10. a is a vertex correction, which has 
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Figure 9: Diagrammatic representation of the collision terms C21 (fig. (a)) and C22 (fig- 
(b)), as computed with the first-order fluctuation 5N^°\k, P). The wavy line with a blob 
denotes the soft gluon propagator in the HTL approximation. The other notations are as 
in Fig. H 




(a) (b) 

Figure 10: Loop corrections to the polarization tensor in Fig. |S|.b, as generated by the 
collision terms depicted in Fig. || 



to supplement the self-energy resummations in Fig. [|b as required by gauge symmetry. 
Similarly, the diagram in Fig. |l0|.b is a different kind of vertex correction which involves 
two hard loops (and one soft one). 

It should be clear by now what are the diagrams generated by further iterations: 
these are the ladder diagrams displayed in Fig. [TT] where the continuous lines are hard 
transverse gluons dressed with the damping rate T, and the (internal) wavy lines are soft 
gluons which may be longitudinal or transverse and are dressed with the hard thermal 
loop. Indeed, the phase-space integrals giving T and Cn are individually dominated by 
soft q momenta, in the range g 2 T ^ q ^ gT. Note also that the diagrams in Fig. |TT] 
may involve an arbitrary number of loops. However, from the previous derivation of 
the Boltzmann equation, we know that, for an external momentum P ~ g 2 T, they all 
contribute at the same order in g, namely, at the same order as the corresponding hard 
thermal loop [|| § (see also Ref. [pfl). 



Except for the colour structure of the vertices, similar diagrams are resummed by the 
Abelian Boltzmann equation as well. There too, the damping rate T of the charged particle 
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Figure 11: A generic ladder diagram of the type resummed in the solution to the Boltz- 
mann equation. 



is dominated by soft momenta transfer, and, in fact, it coincides with the corresponding 
rate for gluons up to a colour factor (compare eqs. ( |2.82| ) and ( |3.107| )). Nevertheless, we 
have seen in Sec. 2.6 that the Abelian collision terms are not dominated by soft momenta 
exchange, but receive contributions from all the momenta between gT and T. In order 
for this to be consistent with the diagrammatic picture in Fig. [TTL there must be some 
cancellations among the Feynman graphs, in such a way that the globl contribution of 
soft (q ^ gT) internal photons cancels out in their sum. 

The relevant cancellations have been discussed in Sec. 2.6, and can be also read off 
the formulae in Sect. 3.8. The contributions of the soft exchanged momenta q ^ gT to 
the individual collision terms are listed in eqs. ( p. 107] ), ( |3.108|) and ( |3.110|) . If there was 
not for the colour structure, the vertex correction in eq. ( |3.108| ) would exactly cancel the 
self-energy correction in eqs. (|3.107|) . This is what happens in the Abelian case (cf. the 
discussion in Sec. 2.6) and explains, for instance, why there is no effect of the damping rate 
7 ~ e 2 Tln(l/e) of the charged particle on the polarization tensor U^ u for soft photons. 
Diagramatically, this corresponds to a cancellation of the soft photon effects in between 
self-energy and vertex corrections in any of the bubbles depicted in Fig. [IT], as it has been 
also verified via the direct analysis of the Feynman graphs, by Lebedev and Smilga [[!(]] 
(see also Refs. g|, g|). 

Similarly, in the absence of colour effects, there would be an exact cancellation 



between the soft contributions to C22 and C23 (cf. eq. (|3.110| )): if we refer to the diagram 
in Fig. []b, this cancellation reflects the well known fact that there is no photon HTL with 
three (or more) external legs fl]. |2|. 

In QCD, on the other hand, these cancellations are not complete, because of the 
non-trivial colour structure (cf. eqs. (|3~T09|) and ( |3~TToD ). In fact, for 5N(k, X) = 5N a T a , 
we have (refer to eq. (|3.108| )) (T a T b ) c d5N c d = (N c /2)6N a b, so that the vertex correction in 
eq. ( |3.108| ) only cancels half of the damping rate contribution in eq. ( |3.107| ). Furthermore, 
since SN = 5N a T a is an antisymmetric colour matrix, the two terms C22 and C23 give 
identical contributions which add together (rather than subtract, as in QED) to give the 
final result in eq. ( |3.110| ). Accordingly, the polarization tensor W^ V {P) for soft (P ~ g 2 T) 
gluons is sensitive to the hard gluon damping rate 7 ~ g 2 T\n(l/g), as discussed at the 
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end of the previous section. (To our knowledge, this has been first observed by Bodecker 

10 

The previous discussion also shows that, for colourless fluctuations in the hot QCD 
plasma (as involved, for instance, in the calculation of the shear viscosity), the pattern 
of the cancellations is the same as in QED: the relevant collision terms are dominated 
by relatively large (gT ^ q ^ T) transferred momenta. Thus, the typical inverse relax- 
ation time for such fluctuations is r _1 ~ g 4 T\n(l/g), where the logarithm comes from 
log(T/m£)), that is, from the screening effects at the scale gT & [14L [15 . 



5 Conclusions 

Starting from first principles, we have derived a Boltzmann equation which describes 
the long wavelength colour excitations of a high temperature Yang-Mills plasma. Our 
derivation relies on a gauge- covariant gradient expansion of the Kadanoff-Baym equations 
for the off-equilibrium dynamics of the plasma. This expansion can be also intrepreted 
as an expansion in powers of g, in the sense that the plasma inhomogeneities (the "soft" 
gradients dx), the strength of the colour mean fields and the off-equilibrium deviations 
of the distribution functions are all contolled by powers of g. Specifically, our derivation 
applies to the case where dx ~ gA ~ g 2 T, so that the soft covariant derivatives Dx = 
dx + igA ~ g 2 T are consistently preserved in the perturbative expansion. This, together 
with a judicious choice of the gauge fixing (the background field gauge) and a proper 
definition of gauge-covariant Wigner functions, has allowed us to maintain explicit gauge 
symmetry at all steps of our construction. 

In this framework, the Boltzmann equation has emerged as the quantum trans- 
port equation at leading order in g. Note that our present, leading order, perturbative 
expansion encompasses several approximations which are generally seen as independent 
approximations when deriving transport equations: the gradient expansion (slowly vary- 
ing disturbances), the weak field expansion (small perturbations), the quasiparticle ap- 
proximation (well defined "quasiparticles" , with a long lifetime) and the (dressed) Born 
approximation (one-gluon exchange) for the collision term. Note also that, in spite of be- 
ing obtained within a systematic perturbative expansion, the Boltzmann equation actually 
resums an infinite series of ladder diagrams and therefore generates a non-perturbative 
effective theory for the "ultrasoft" (dx ~ g 2 T) colour fields. For instance, the polarization 
tensor for these fields, as obtained by solving the Boltzmann equation, is equivalent to 
the sum of an infinite number of ladder diagrams of the perturbation theory. 

Remarkably, the colour structure of the collision terms is precisely that predicted 
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by Arnold, Son and Yaffe, on the basis of simple heuristic arguments |7| (this is also 



consistent with some previous results in Refs. [17], [10], [12|])- Accordingly, the effective 
theory generated by our Boltzmann equation is precisely that obtained by Bodecker using 
a different method |J (see Ref. for a discussion of the relation between the Boltzmann 
equation and Bodecker's theory). 

It is important to emphasize the accuracy limits of the Boltzmann equation pre- 
sented here. The collision term in eq. ( |1.2| ) is only known to logarithmic accuracy, which 
means that both its functional form and the overall coefficient 7 which fixes the time scale 
for colour relaxation are known only up to corrections of relative order l/ln(l/g). For 
instance, 7 = aN c T(\n(l/ g) + 0(1)) where the constant terms under the logarithm is not 
accurately given by the present formalism. There are two main sources for such a limita- 
tion: (i) The gradient expansion in the presence of long range interactions: the one-gluon 
exchange interaction has a typical range ~ 1/q with g 2 T ^ q ^ gT which is marginally 
the same as the inhomogeneity scale in the problem, A ~ l/g 2 T. Since the collision terms 
are only logarithmically sensitive to the low momenta q ~ g 2 T, it follows that the gradient 
expansion for the collision kernel (i.e., for the exchanged gluon propagator in eq. fl3.90| )) 
is only correct to logarithmic accuracy. As discussed after eq. (|3.92|) , this limitation can 
be avoided, at least in principle, by relaxing the gradient expansion so as to allow for a 
collision kernel which is non-local in X. (ii) Still related to the long-ranged collision ker- 
nel: the behaviour of the gluon propagator at very soft momenta q ~ g 2 T is not correctly 
described by the HTL approximation, since it is sensitive to the non-perturbative physics. 

Of course, the Boltzmann equation constructed in this paper can also be used to 
study colourless fluctuations, as involved, e.g., in the calculation of the shear viscosity [§, 
|14], j34j. The relevant collision terms are displayed in eqs. ( |3.99|) , (|3.101|) and (|3.104j) , where 
the off-equilibrium fluctuations have now a trivial colour structure (e.g., 5G^ d (p, X) = 
3dd$G(p, X), etc.). This yields the standard Boltzmann equation (in its linearized version), 
as already used in calculations of the shear viscosity for the quark gluon plasma || |14}| . 
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